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ABSTBACT 


Si 

Pointing  error  standard  deviations  for  two  theodolites, 
the  Wild  T-2  and  Odoa  Aztrac,  were  determined  under  condi¬ 
tions  closely  approximating  those  of  range-azimuth  or 
azieuth-aziauth  hydrographic  surveys.  Pointing  errors  found 
for  both  instruments  were  about  1.3  meters,  and  were  inde¬ 
pendent  of  distance.  No  statistical  difference  between  the 
errors  of  the  two  instruments  was  found.  The  accuracy  of 
the  interpolation  methods  used  by  the  National  Ocean  Service 
(NOS)  for  range-azimuth  positioning  were  investigated,  and 
an  average  inverse  distance  of  about  2.5  meters  was  observed 
between  interpolated  positions  and  corresponding  observed 
positions.  The  overall  range-azimuth  position  errors  of  the 
two  theodolites  were  then  compared  to  positioning  standards 
of  NOS  and  the  International  Hydrographic  Organization, 
using  assumed  ranging  standard  deviations  of  1.3  and  3.0 
meters.  Both  instruments  met  all  standards  except  the  NOS 
range-azimuth  standard  for  1:5,000  scale  surveys. 
Interpolated  positions  may  fail  to  meat  more  of  the  stan¬ 
dards  because  of  additional  inherent  error. 


T1BLE  OP  CONTENTS 


I.  INTRODUCTION . 9 

A.  THE  RANGE-AZIMUTH  POSITIONING  METHOD . 9 

B.  HYDROGRAPHIC  POSITION  ERROR  STANDARDS  ....  12 

C.  OBJECTIVES . 14 

II.  ERROR  INDICES  AND  RANGE-AZIMUTH  GEOMETRY  .  20 

A.  DEFINITIONS . 20 

1.  Slanders . 22 

2.  Systematic  Errors . 23 

3.  Random  Errors  . . 24 

3.  TWO -DIHEN SION AL  ERROR  FIGURES  .  26 

1.  Concentric  and  Eccentric  Geometry  ....  27 

2.  The  Error  Ellipse . 29 

3.  Root  Mean  Square  Distance . 33 

4.  circular  Standard  Error . 35 

C.  THE  ERROR  OF  AN  INTERPOLATED  FIX . 37 

1.  Interpolation  Algorithms . 37 

2.  Error  Propagation  .  39 

III.  EXPERIMENT  DESIGN  AND  IMPLEMENTATION  .  43 

A.  FIELD  WORK . 43 

B.  THEODOLITE  POINTING  ERROR  .  47 

C.  INTERPOLATION  ALGORITHM  EVALUATION . 49 

D.  CHOICE  OF  EXPERIMENTAL  CONDITIONS  .  51 

IV.  RESULTS  AND  DATA  ANALYSIS . 56 

A.  DATA  PROCESSING  SYSTEM . 56 

B.  POINTING  ERROR  DETERMINATION  .  56 

C.  INTERPOLATION  EVALUATION . 64 


D.  ANALYSIS  OF  FACTORS  AFFECTING  THE  RESULTS  .  .  66 


5 


B.  APPLICATION  TO  POSITION  ERROR  STANDARDS  .  .  . 


V.  CONCLUSIONS  AND  RECOMMENDATIONS . 79 

APPENDIX  A:  MEANS  AND  STANDARD  DEVIATIONS  OP 

ACQUIRED  DATA  SETS . 85 

APPENDIX  B:  3EODETIC  POSITION  OP  HORIZONTAL  CONTROL 

STATIONS . 91 

APPENDIX  C:  EMPIRICAL  PROBABILITY  DENSITY  FUNCTION 

PLOTS . 92 

BIBLIOGRAPHY  .  98 

INITIAL  DISTRIBUTION  LIST . 101 


6 


LIST  OP  TABLES 


Circular  Error  Formulae . 36 

Data  Acquisition  Sequence  of  Brants . 45 

Pointing  Error  standard  Deviation  (pooled 

estimates)  .  ....58 

Summary  of  AHOVA  Results  at  95*  Confidence  ....  63 

Results  of  Interpolation  Evaluation  .  65 

Corrected  Aztrac  Standard  Deviation . 70 

Experiment  Precision  and  Theodolite  Error 

Position  Standards  Comparison . 

3000  deter  Arc  . . 

1500  deter  Arc . 

1000  deter  Arc . 

700  deter  Arc  . 

500  Heter  Arc  . 

300  det9r  Arc . 


LIST  OP  FIGOBES 


1.1  Illustration  of  Range-  Aziiuth  Positioning  ...  10 

1.2  Illustration  of  Aztrac  Shore  Station . 17 

1.3  The  Aztrac  Transmitting  Unit . 19 

2.1  The  Horsal  Probability  Curve . 25 

2.2  Eccentric  Range-Azimuth  Geometry . 27 

2.3  Concentric  Sange-A ziauth  Seoaetry . 28 

2.4  Eccentric  Error  compensation  .  29 

2.5  Error  Ellipse  and  dr«,* . 31 

2.6  A  Range-Azimuth  Position . 32 

2.7  Variation  of  drw,s  Probability . 34 

2.8  Example  of  Angular  Interpolation  .  38 

3.1  Sketch  of  the  survey  Area . 44 

3.2  Lines  of  Position  Observed  to  the  vessel  ....  45 

3.3  Oncertainty  of  an  observed  Error . 48 

3.4  Interpolation  of  Angle  Only . 50 

3.5  Interpolation  of  Angle  and  Distance  .  50 

4.1  One -Dimensional  Difference  Between  Means  ....  66 

4.2  Two-Dimensional  Difference  Between  Weans  ....  67 

4.3  Probability  Density  vs.  Error . 68 

4.4  Results  Compared  to  dr«s . 74 

4.5  Results  Compared  to  90*  Probability . 76 

C.  1  Probability  Density  Plot:  300  m  arc  . . 92 

C.2  Probability  Density  Plot:  500  a  arc . 93 

C.3  Probability  Density  Plot:  700  a  arc . 94 

C.4  Probability  Density  Plot:  1000  a  arc . 95 

C.5  Probability  Density  Plot:  1500  a  arc . 96 

C.6  Probability  Density  Plot:  3000  a  arc . 97 


8 


I.  IXH&B2SZI21 


k.  THE  HiHGB-iZIHUTH  P0SIT20VIHG  BBT  HOD 

The  fundaaental  purpose  of  a  hydrographic  survey  is 
defined  by  Inghaa  (1974)  as  being  to  "depict  the  relief  of 
the  seabed,  including  all  features,  natural  and  aanaade,  and 
to  indicate  the  nature  of  the  seabed  in  a  Banner  sinilar  to 
the  topographic  sap  of  land  areas."  He  goes  on  to  describe 
two  factors  defining  a  single  point  on  the  seabed: 

(i)  "The  position  of  the  point  in  the 
horizontal  plane  in,  for  exaaple,  latitude 
and  longitude,  grid  co-ordinates  or  anal es 
and  distances  froa  known  control  points. 


(ii)  The  depth  of  the  point  below  the  sea 
surface,  corrected  for  the  vertical 
distance  between  the  point  of  aeasurenent 
and  water  level  and  for  the  height  of  the 
tide  above  the  datus  or  reference  level  to 
which  depths  are  to  be  related." 

Thus  the  hydrographer  aust  answer  the  two  priaary  ques¬ 
tions  of  "how  deep"  and  "where"  for  each  of  the  thousands  of 
soundings  acquired  on  every  survey.  Because  every  area  to 
be  surveyed  has  different  geophysical  characteristics  and 
levels  of  use,  the  hydrographer  aust  possess  a  suite  of 
tools  and  techniques  to  accoaplish  aach  survey,  A  survey  of 
a  large  aetropolitan  harbor  requires  different  equipaent  and 
aeasureaent  precision  than  one  for  a  deep  ocean  area. 

only  the  first  of  Inghaa' s  two  factors  cited  above  is 
considered,  and  it  is  further  narrowed  in  scope  to  techni¬ 
ques  used  in  the  nost  precise  surveys.  Such  a  survey  sight 
be  one  of  a  winding,  narrow  river  carrying  deep  draft 
vessels,  or  perhaps  a  very  large  scale  survey  of  an  inner 
harbor.  Both  areas  require  the  highest  positioning  accuracy 
and  a  ainiaua  of  shore  control  stations. 
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Any  ■•thod  of  positioning  esploys  the  intersection  of 
lines  of  position  (LOP's)  to  construct  a  fix.  Although 
advanced  eethods  say  use  sultiple  LDP's,  traditional  hydrog¬ 
raphy  uses  the  simple  intersection  of  two  lines  to  fix  the 
vessel's  position.  The  vessel  is  located  somewhere  along 
each  of  two  lines  of  position#  and  the  only  point  satisfying 
these  conditions  is  the  intersection  of  the  lines. 

The  error  associated  with  one  of  the  simplest  posi¬ 
tioning  methods,  that  of  the  two  LDP  range-azimuth  fix# 
which  is  illustrated  in  Figure  1.1,  is  analyzed. 


Also  called  the  rho/theta  method#  range-azimuth  positioning 
consists  of  the  observation  of  a  distance  and  an  azimuth  to 
a  vessel  from  either  one  or  two  known  locations  (Ombach# 

1976  ].  An  example  of  this  method  is  the  use  of  radar  aboard 
ship.  A  relative  position  for  a  radar  contact  is  determined 
by  observing  a  radar  range  and  azimuth#  or  a  radar  range  and 
visual  azimuth#  to  a  contact.  The  two  lines  of  position 
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always  intersect  at  right  angles  because  the  observation  is 
■ad*  fro*  a  single  point,  and  this  concentric  geometry 
provides  the  strongest  fix  possible.  Mariners  also  know 
that  the  fix  obtained  via  a  visual  aziauth  is  stronger  than 
the  one  using  a  radar  azimuth,  because  the  visual  bearing  is 
■ore  accurate. 

This  example  shows  the  advantages  that  Bake  the  range- 
aziauth  method  popular  for  hydrography.  It  provides  the 
geometrically  strongest  possible  fix,  and  only  one  location 
on  shore  need  be  occupied  to  control  the  survey.  Such  a 
positioning  method  is  ideal  in  harbors  or  rivers  where 
aaxiaua  accuracy  is  needed  but  where  obstructions  make  ether 
types  of  fix  geoaetry  impractical.  In  1982  the  o.S. 

National  Ocean  Service  (NOS)  obtained  twenty  thousand  linear 
nautical  miles  of  launch  hydrography,  and  sixty  percent  of 
this  was  controlled  by  the  range-azimuth  method  [Wallace, 
1983  ]. 

There  are  limitations  associated  with  this  method  just 
as  with  any  fix  geoaetry.  It  is  labor  intensive  and 
requires  more  radio  communication  (to  establish  fix  timing) 
than  mest  other  methods.  In  totally  nonautomated  situ¬ 
ations,  distances  to  the  survey  vessel  are  recorded  manually 
aboard  the  vessel,  and  azimuths  are  recorded  ashore  by  the 
theodolite  observer  at  prescribed  intervals.  These  fix  data 
are  later  put  into  computer  compatible  digital  form  via  a 
process  called  logging.  Manual  recording  of  these  fix  data 
are  generally  toe  slow  to  position  every  sounding. 

Therefore,  individual  sounding  positions  must  be 
interpolated  from  the  observed  fixes. 

Systeas  have  been  designed  that  have  an  intermediate 
level  of  automation.  The  NOS  Hydroplot  System  is  an  example 
of  this  type  'Wallace,  1967  ].  when  used  in  the  range- 
aziauth  aode,  the  vessel  is  usually  steered  along  arcs  of 
constant  range  froa  the  theodolite  station,  and  the 
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Hydroplot  System  automatically  records  a  distance  measure¬ 
ment  for  each  sounding.  Aziauths  are  not  teleaetered  to  the 
vessel  but  are  relayed  over  voice  radio  and  are  manually 
entered  into  the  computer  system.  Since  the  maximum  data 
rate  is  about  two  angles  per  minute,  the  interpolation  of 
angles  for  sounding  positions  between  fixes  is  necessary. 

Recently  a  digital  theodolite,  the  Odom  Aztrac,  has  been 
developed  which  can  record  and  telemeter  angles  with  great 
speed  —  up  to  ten  angles  per  second  [Odom  Offshore  Surveys, 
Inc.,  1982].  A  computer  system  aboard  the  survey  vessel  can 
thus  record  and  plot  an  observed  position  for  each  sounding. 
This  rapid  position  fixing,  combined  with  a  computer's 
ability  to  provide  cross-track  error  indications  to  the 
helmsman,  enables  the  hydro grapher  to  systematically  cover  a 
survey  area  with  maximum  efficiency  by  running  straight  and 
parallel  sounding  lines. 

The  Aztrac  system  is  still  considered  a  seeiautomated 
system  because  an  observer  is  required  to  manually  track  the 
vessel  with  the  theodolite.  Two  fully  automated  range- 
azimuth  systems  which  feature  fully  automatic  tracking  have 
been  devolcped.  One  is  the  Polarfix  system  developed  by 
Krupp-Atlas  Elektronik  in  Germany  [Smith,  1983],  and  the 
other  is  the  Artemis  system  developed  by  Christiaan 
Huygenslaborator ium  in  the  Netherlands  [Newell,  1981]. 

B.  HYDROGRAPHIC  POSITION  ERROR  STANDARDS 

Historically,  most  national  hydrographic  organizations, 
as  well  as  the  International  Hydrographic  Organization 
(IHO) ,  have  used  linear  plotting  error  at  the  scale  of  the 
survey  to  be  the  standard  for  sounding  position  accuracy. 
Prior  to  1982,  the  standards  recommended  by  IHO  [IHO,  1968] 
were: 


"The  indicated  repeatability  of  a  fix 
(accuracy  of  location  referred  to  shore 
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control)  in  the  operating  area,  whether 
observed  by  visual  or  electronic  methods, 
ccsbined  with  plotting  error,  shall  seldom 
exceed  1.5  mm  (0.05  in)  at  the  scale  of 
the  survey." 


The  I HO  recently  published  new  recommendations  for  error 
standards  [IHO,  1982]  which  are: 

"...  any  probable  error,  measured  relative 
tc  shore  control,  shall  seldom  exceed 
twice  the  minimum  plottable  error  at  the 
scale  of  the  survey  (normally  1.0  mm  on 
paper)  . " 

Neither  of  the  IHO  standards  make  any  reference  as  to 
what  probability  level  they  apply.  Hunson  (1977)  inter¬ 
preted  the  words  "shall  seldom  exceed"  in  the  above  state- 
ments  to  mean  "less  than  10 %  of  the  time",  which  seems 
reasonable.  The  1982  IHO  standard  is  somewhat  confusing  due 
to  its  use  cf  both  the  terms  "seldom  exceed"  and  "probable 
error".  The  latter  tern  is  associated  with  a  50X  prob¬ 
ability  by  most  statisticians  including  Green  wait  (1971). 
However,  the  author  of  these  standards.  Commodore  A.  H. 
Cooper,  HAN  (retd)  ,  has  stated  that  he  intended  no  statis¬ 
tical  significance  to  the  term  "probable  error"  [Wallace, 
1983 ]. 

The  NOS  has  not  yet  incorporated  the  latest  IHO  stan¬ 
dards,  but  such  action  is  being  considered  in  some  fora 
(Wallace,  1983  ].  Current  NOS  standards  have  been  developed 
to  ensure  that  "accuracies  attained  for  all  hydrographic 
surveys  conducted  by  NOS  shall  equal  or  exceed  the  specifi¬ 
cations"  of  the  1968  IHO  standards  [Uabach,  1976].  Unlike 
the  international  standards,  the  NOS  standards  for  all  elec¬ 
tronic  positioning  systems  use  the  concept  of  root  mean 
square  error  (dr«s  or  rmse)  ,  which  has  a  somewhat  variable 
probability  of  between  6  8.3  and  63.2  percent.  The  NOS  stan¬ 
dards  for  fully  visual  and  for  hybrid  (combination  elec¬ 
tronic  and  visual)  positioning  have  no  explicit  reference  to 
probability. 


Specific  operational  standards  for  range- aziauth  posi¬ 
tioning  have  been  neglected  by  aany  hydrographic  organiza¬ 
tions.  However,  nos  [(Jabach,  1976]  requires  the  following 
observational  procedures  be  followed  for  all  range-azieuth 
positions. 

"Objects  sighted  on  should  be  at  least  500 
■  from  the  theodolite...  the  aziauth 
check  should  not  exceed  one  ainute  of 
arc...  observed  aziauths  or  directions  to 
the  sounding  vessel  for  a  position  fix 
shall  be  read  to  the  nearest  1  ain  of  arc 
or  better  if  necessary  to  produce  a  posi¬ 
tional  accuracy  of  0.5  aa  at  the  scale  of 
the  survey." 

Since  the  range-aziauth  aethod  is  classified  as  a  hybrid 
positioning  systea,  it  is  not  referenced  tc  any  particular 
probability,  but  a  reasonable  assuaption  aay  be  aade  that 
the  drw«  concept  also  applies  in  this  case. 

The  U.S.  Naval  Oceanographic  Office  (NAVOCEANO)  also 
requires  that  its  surveys  aeet  the  standards  of  the  nos 
Hydrographic  Hanual.  The  Aray  Corps  of  Engineers  presently 
have  no  formal  positioning  requirements  that  aust  be  aet  by 
all  districts,  although  draft  specifications  are  being 
written  at  this  tiae  [Hart, 1983],  The  range-aziauth  tech¬ 
nique  and  its  applicability  to  Corps  of  Engineers  surveys  is 
discussed  in  Hart  (1977).  No  specific  reguireaents  for 
range-aziauth  positioning  could  be  found  for  either  the 
Canadian  Hydrographic  Service  or  ths  British  Hydrographic 
Service.  Palikaris  (1983)  also  reports  no  published  stan¬ 
dards  for  these  organizations. 

C.  OBJECTIVES 

All  position  error  standards  using  an  explicit  prob¬ 
ability  are  based  on  the  idea  that  an  observation  is  a 
normally  distributed  random  variable  with  zero  Bean  and 
standard  deviation  <r  .  These  standards  require  a  value  for 
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the  standard  deviation  of  the  conponent  lines  of  position 
that  sake  up  the  fix.  The  standard  deviation  is  a  value 
such  that  there  is  a  68.27*  probability  of  an  observation 
falling  within  ±1  <T  of  the  aean.  It  is  unfortunate  that 
often  the  1  <T  values  of  a  hydrographic  measurement  are 
simply  not  known,  or  known  only  for  ideal  conditions  and 
provided  by  manufacturers  who  have  a  vested  interest  in  the 
measuring  instrument. 

This  paper  and  experiment,  then,  has  as  its  primary 
purpose  the  determination  of  a  pointing  error  standard  devi¬ 
ation  for  two  theodolites  used  to  measure  azimuths  under 
hydrographic  conditions.  One  instrument  (the  Wild  T-2)  is 
the  standard  used  by  NOS  field  units.  The  other  (the  Odom 
iztrac)  is  a  new  digital  talematering  theodolite,  which  may 
prove  useful  in  automating  the  presantly  tedious  manual 
methods.  No  astimation  of  standard  deviation  has  ever  been 
made  on  these  devices  under  typical  range-azimuth 
conditions. 

The  T-2  is  the  standard  instrument  used  by  nos  surveying 
parties  for  land  surveying  and,  to  a  lesser  extent,  hydrog¬ 
raphy.  It  is  a  very  precise  instrumant  used  in  third-order 
horizontal  control  and  provides  an  angular  resolution  of  one 
second  of  arc.  However,  it  dees  possess  features  that  are 
lass  than  ideal  for  range-azimuth  work.  The  horizontal 
tangent  screw  of  this  instrument  is  awkward  for  range- 
azimuth  surveying  because  it  is  not  an  infinite  gearing 
device.  The  observer  often  encountars  the  end  of  the  drive 
mechanism,  stopping  the  instrument's  movement  while  tracking 
the  vessel,  k  solution  to  thi3  problen  for  many  observers 
is  to  track  the  vessel  with  the  tangent  screw  unclasped, 
then  clamp  the  screw  only  seconds  before  the  fix  occurs. 

The  inverted  image  feature  of  the  T-2's  in  use  by  NOS 
simplifies  the  optical  system  and  reduces  optical  error. 

This  is  satisfactory  for  land  survey  work  but  creates  some 


confusion  when  tracking  a  fast  soring  vessel.  This  is 
because  the  observer  sees  an  iaage  of  the  survey  boat  upside 
down  and  soving  apparently  in  the  opposite  direction  from 
its  actual  aoveaent.  Another  disadvantage  is  that  the  mech- 
anisa  for  reading  the  horizontal  circle  of  a  T-2  is  sore 
coaplicated  than  desired  for  the  rapid  observations  neces¬ 
sary  in  hydrographic  survey  work.  Ihe  operator  is  required 
to  step  tracking  the  vessel ,  reaove  his  eye  froa  the  tele¬ 
scope,  and  use  an  auxiliary  eyepiece  to  read  the  angle. 

There  is  a  practical  liait  of  about  30  seconds  to  the  speed 
with  which  successive  angles  can  be  observed  and  read.  The 
particular  unit  tested  was  serial  number  30504. 

The  Odoa  Aztrac  theodolite  is  a  seal-automated, 
line-of-sight  angle  aeasuring  systaa.  The  Aztrac  system 
consists  of  a  wild  T-16  theodolite  (serial  number  2534880 
was  tested)  which  was  modified  for  infinite  tangent  drive 
and  to  provide  angular  information  in  a  digital  format.  The 
shore  unit  decodes  the  observed  angle,  determines  the  direc¬ 
tion  of  rotation  of  the  instrument  and  displays  the  angle  or. 
it3  front  panel.  The  angle  is  then  converted  to  binary 
coded  decimal  (BCD)  format  and  used  to  frequency  shift  key 
(PSK)  an  PM  transmitter  to  link  tha  data  with  the  survey 
vessel.  The  data  is  transmitted  at  the  rate  of  10  angles 
per  second.  On  board  the  vessel  the  Aztrac  receiver 
converts  the  received  data  to  parallel  form  and  displays  it 
on  the  front  panel  for  nanual  recording.  Pot  automated 
recording  or  processing  by  onboard  computer  a  serial  data 
ouout  is  provided.  The  Aztrac  has  an  angular  resolution  of 
0.01  degree  (36  arc  seconds)  (Odoa,  1983]. 

With  the  notable  exception  of  angular  resolution,  this 
theodolite  is  sore  appropriate  for  range-aziauth  hydrography 
than  the  wild  T-2.  It  has  an  erect  image  and  infinite 
tangent  screw  which  allow  the  vessel  to  be -constantly 
tracked.  Its  digital  output  requires  no  action  on  the 
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operator's  part  to  record  observed  angles.  An  additional 
positive  feature  of  the  digital  output  is  the  ability  to 
rapidly  zero  the  instrunent  on  an  initial  pointing.  This  is 
done  by  pressing  a  aanual  reset  button  on  the  instrument 
control  panel.  Figures  1.2  and  1.3  show  photographs  of  the 
Aztrac  equipaent.  Figure  1.2  illustrates  a  typical  shore 
station.  The  observer  is  adjusting  the  Aztrac  instrument, 
and  the  transaitter  unit  is  on  the  ground  to  the  right  of 
the  Aztrac  tripod.  The  distance  measuring  equipaent  is 
mounted  on  a  tripod  behind  the  observer,  and  the  Aztrac 
transmitting  antenna  is  at  the  top  of  a  pole  on  the  extreme 
right  of  the  photograph.  Pigure  1.3  shoes  the  Aztrac  trans¬ 
aitter  unit,  with  its  digital  angular  display  in  hundredths 
of  a  degree.  Both  photographs  were  provided  by  Odom 
Offshore  Surveys,  Inc. 

A  seccnd  objective  of  this  paper  is  to  evaluate  the 
interpolation  methods  used  by  the  NDS  for  range-azimuth 
work.  The  availability  of  the  digital  theodolite,  with  its 
direct  measurement  of  all  positions,  enabled  a  comparison  of 
observed  and  interpolated  fixes  to  be  made.  This  made 
possible  an  estimate  of  whether  interpolated  positions  meet 
required  accuracy  standards,  and  whether  direct  measurement 
of  all  positions  is  needed.  The  estimate  presented  here  was 
not  made  statistically  rigorous  so  that  the  thesis  could  be 
kept  to  a  manageable  size.  More  theoretical  statistical 
work  is  needed  to  fully  reduce  the  interpolation  data. 

The  final  objective  of  this  investigation  is  to  compare 
the  position  errors  of  these  two  instruments  with  the 
various  position  error  standards  discussed  in  section  B  of 
this  chapter.  The  conclusions  resulting  from  this  objective 
will  assist  the  hydrogr&pher  to  select  equipaent  and  oper¬ 
ating  conditions  that  seat  required  position  error 
standards. 
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XI.  BUM  IflJ2X£ZS  Ufi  SUSfciUASXi  SSSflSXfiJ 

The  method  of  range-aziauth  positioning  is  asaally 
selected  for  large  scale  surveys  because  of  its  siaplicity 
and  accuracy.  This  is  the  result  of  both  angle  and  distance 
aeasureaent  devices  being  co-locatei,  with  the  intersection 
of  the  two  lines  of  position  always  being  ninety  degrees. 

In  practice,  however,  the  co-location  of  both  instruaents  is 
often  not  achieved.  The  result  is  an  eccentric  geoaetry  for 
the  position  fix.  This  section  will  analyze  the  geoaetry  of 
both  eccentric  and  concentric  fixes.  Position  error  indices 
in  coaaon  use  will  be  reviewed  and  analyzed  for  the  special 
cases  of  range-aziauth  methods,  and  the  error  of  an 
interpolated  fix  will  be  derived. 

1.  DBFI1ITI0HS 

Although  a  complete  and  general  treataent  of  error 
theory  will  not  be  presented,  soae  basic  definitions  are 
necessary  to  understand  the  data  analysis  presented  here. 

The  ideas  in  this  section  are  included  in  aany  basic  statis¬ 
tics  textbooks,  and  were  specifically  drawn  froa  Wonnacott 
(1977),  Bowditch  (1977),  Heinzen  (1977),  Kaplan  (1980),  and 
Davis  (1981). 

Error  aay  be  defined  as  "the  difference  between  a 
specific  value  and  the  correct  or  standard  value"  [Bowditch, 
1977],  or  as  "the  difference  between  a  given  aeasureaent  and 
the  "true"  or  "exact"  value  of  the  measured  quantity" 

[Davis,  1981].  Batheaatica lly  it  can  be  defined  as: 


where  •  is  the  error,  xx  is  an  observation,  and  T  is  the 
"correct"  or  "true"  value.  The  word  error  implies  that 
thars  is  a  known  true  Talus  for  a  quantity,  with  which  a 
ssasursssnt  say  b#  compared  to  find  tha  "error"  associated 
with  that  aaasursaant.  Sines  tha  true  talus  of  a  measured 
quantity  is  rarely  knewn,  ths  tern  "error"  is  not  precisely 
correct.  Davis  (1981)  states  that  it  is  aore  appropriate  to 
speak  of  ths  theory  of  observations  rather  than  the  theory 
of  errors,  but  it  can  be  shown  that  the  difference  between 
ths  two  is  largely  one  of  sesantics. 

A  single  isasuresent  of  a  particular  quantity  aay  be 
considered  sufficient  for  aany  purposes,  even  if  it  is  known 
that  additional  aeasureaents  will  probably  be  slightly 
different  than  the  first.  If  ths  quantity  to  be  measured  is 
of  sufficient  inport ance,  then  aultiple  aeasureaents  are 
aads  and  the  saaple  aean,  Z,  is  used.  Each  of  these 
aultiple  aeasureaents  can  be  a  considered  numerical  value 
for  a  random  variable.  A  random  variable  is  one  that  takas 
on  a  range  of  possible  values,  each  associated  with  a 
particular  probablility. 

The  saapla  aean  aay  be  expressed  aatheaat ically  by  egua- 
tion  2.2  [Wcnnacott,  1977]: 

—  i 

^  *  n  (2.2) 

where  n  is  the  saaple  size.  If  the  saaple  size  were 
increased  without  limit  (n  —*•<»),  aquation  2.2  would  give 
the  population  aean .  The  saaple  aean  is  always  an  esti¬ 
mate  of  the  population  aean,  which  is  never  directly 
computed.  This  leads  to  the  concept  of  the  residual,  v, 
which  is  the  difference  between  the  estiaate  7  of  the  popu¬ 
lation  aean  and  the  observation  x^  .  This  is  shown  in 
equation  2.3. 


(2.3) 


y  *  X  - 

The  residual  is  cosputa tionally  the  negative  of  the 
error.  Nevertheless,  equation  2.3  is  sore  appropriate 
because  it  uses  an  estimate,  7,  of  the  unknowable  population 
aean .  The  presence  of  X  in  equations  2.3  iaplies  that 
multiple  eeasureaents  have  been  sale,  and  allows  a  partic¬ 
ular  confidence  to  be  assigned  to  the  estiaate  of 
depending  on  the  nuiber  of  such  seasuresents.  Because  the 
word  error  is  still  used  in  such  of  the  hydrographic  profes¬ 
sion,  it  will  be  used  interchangeably  in  this  paper  with  the 
ten  residual.  It  is  iaportant,  however,  to  understand  that 
the  concept  of  the  residual,  whatever  its  nase  nay  be,  is 
fundasental  to  any  aeasureaent  operation. 

Errors  are  classically  divided  into  three  groups:  blun¬ 
ders,  systematic  error,  and  randos  errors  [Greenwalt,  1962]. 
Eowditch  (1977)  and  Davis  (1981)  do  not  classify  errors  as 
including  blunders,  but  like  the  tarn  error  itself,  the 
distinction  is  largely  a  semantic  one.  Ideally  blunders  and 
systenatic  errors  are  completely  eliminated  from  the  data. 
The  most  precise  measurements  reduce  random  error  as  much  as 
possible,  but  it  can  never  be  completely  eliminated. 

i .  iianiats 

Blunders  are  mistakes  that  are  "usually  gross  in 
magnitude  compared  to  the  other  two  types  of  errors"  [Davis, 
1981],  and  are  most  often  caused  by  carelessness  on  the  part 
of  the  observer,  or  by  grossly  malfunctioning  observing 
equipment.  They  are  usually  detected  and  eliminated  by 
procedural  checks  during  the  data  asquistion  process.  The 
recognition  of  a  blunder  is  not  always  easy,  since  a  blunder 
"may  have  any  magnitude,  and  nay  be  positive  or  negative" 
(Bowditch,  1977]. 
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Systeaatic  errors  are  defined  by  Davis  (1981)  as 
those  that  ocear  "according  to  a  SY5t9i  which*  if  known*  can 
always  be  expressed  by  sathesatical  forsulation. "  This 
aathesatical  aodel  results  in  correctors  that  are  applied  to 
all  seasureaents  obtained*  thus  eliainating  the  systeaatic 
errors  froa  the  observations.  The  aodel  aay  be  as  siaple  as 
a  constant  corrector  subtracted  froa  lengths  obtained  with  a 
steel  tape*  or  it  aay  be  as  coaplicated  as  aoaelling  the 
effects  of  ataospheric  refraction  on  electronic  distance 
aeasuring  equipaent. 

If  the  systeaatic  error  is  such  that  it  cannot  be 
aodelled*  it  is  then  estiaated  by  a  process  known  as  cali¬ 
bration.  Kaplan  (1980)  defines  calibration  as  the  process 
of  coaparing  the  aeasuring  instruaent  against  a  "known" 
standard.  The  word  "known"  is  usually  operationally  defined 
as  a  aeasureaent  operation  or  instruaent  that  is  much  aore 
accurate  than  the  one  being  calibrated.  The  difference 
between  the  observed  and  standard  value  is  used  as  an  esti- 
aate  of  the  total  effect  of  all  systeaatic  errors  present. 
This  process  is  very  clcse  to  the  classical  concept  of 
"errors"  presented  above*  and  is  entirely  proper  for  use  in 
the  correction  of  systeaatic  errors  [Davis*  1981].  of 
course*  one  aust  be  careful  to  apply  the  corrector  only  to 
those  aeasuraents  aade  under  the  saae  conditions  as  the 
calibration. 

k  systeaatic  error  found  in  theodolite  or  sextant 
observations  is  known  as  the  personal  error  of  the  observer 
[Hueller*  1969],  [Bowditch*  1977].  This  type  of  error  is 
rarely  quantified  for  hydro gaphic  applications*  but  never¬ 
theless  it  does  exist.  The  observer  aust  rely  on  the  senses 
of  hearing  and  vision  to  sake  aeasureaents*  which  vary 
between  individuals  as  well  as  with  tiae  in  one  individual. 


some  personal  errors  are  constant  and  soae  are  erratic 
[Davis,  1981].  These  errors  are  ainiaized  by  training  and 
standardizing  observational  procedures.  The  best  way  to 
eliainate  personal  error  is  by  the  use  of  coapletely 
autoaated  observation  equipment. 

3.  BiaiSl  LH2I S 

"Fandom  errors  are  chance  errors,  unpredictable  in 
aagnitude  or  sign",  and  are  "governed  by  the  lavs  of  prob- 
ability"  [Bovditch,  1977  ].  If  one  assumes  that  all  blurders 
and  systeaatic  errors  have  been  removed  from  the  observa¬ 
tions,  the  remaining  values  can  be  regarded  as  sample  values 
for  a  random  variable.  As  noted  earlier,  a  random  variable 
can  take  on  a  range  of  values,  each  associated  with  a 
particular  probability.  A  random  arror  has  high  probability 
of  being  close  to  the  population  aean,^  ,  and  a  low  prob¬ 
ability  of  being  very  much  different  than ^ m  [Greenvalt, 

1962  ]. 

A  probability  density  function  expresses  the  rela¬ 
tion  between  a  value  for  a  random  variable  and  the  prob¬ 
ability  of  its  occurrence.  Hydrographic  survey  measurements 
often  use  the  normal  or  Gaussian  probability  density  func¬ 
tion.  A  concise  explanation  of  this  function  is  given  in 
Greenwalt  (1962)  and  Kaplan  (1980).  The  function  itself  is 
given  as  equation  2.4,  where  p  (v)  is  the  probability  of  the 
occurrence  of  a  particular  residual  v,  and  <T  is  the  popula¬ 
tion  variance  which  is  approximated  by  the  sample  variance, 
S\  given  by  equation  2.5. 

* 

(2.5) 


2<r' 


(2.4) 
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To  find  the  probability  of  a  residual  falling 
between  two  residuals  v,  and  v4  ,  equation  2.4  oust  be  inte¬ 
grated  over  the  interval  v,  to  v* .  This  corresponds  to  the 
area  under  the  gaussian  curve  between  those  two  points,  as 
is  shewn  in  Figure  2.1.  If  p(v)  ware  integrated  froe  -1<r  to 
♦  1  <r,  the  area  under  the  curve  would  be  68.27%  of  the  total 
area.  This  aeans  that  there  is  a  68.  27%  probability  of  a 
particular  residual  falling  between  plus  or  minus  one  stan¬ 
dard  deviation  of  the  lean,  where  the  standard  deviation,  <T , 
is  defined  as  the  square  root  of  the  variance  given  in 
equation  2.  5. 


ply) 


Figure  2.1  The  Sormal  Probability  Curve. 

The  combined  effect  of  blunders,  systematic  errors 
and  randem  errors  can  now  be  seen  in  overview.  Xf  it  is 
assumed  that  blunders  and  systematic  errors  have  been 
completely  eliminated  from  a  set  of  observations,  there 
remain  only  random  errors.  If  the  sample  size  is  large 
enough,  then  the  sample  mean  and  variance  are  good 


approximations  of  the  population  naan  and  variance.  That 
is,  there  is  a  68.27*  probability  that  any  future  measure- 
■ents  aade  under  the  saae  conditions  will  not  fall  farther 
than  plus  or  ainus  one  standard  deviation  from  the  mean. 

B.  T  VO -DIHENS  IONA  L  BBBOR  PIGOBES 

The  two-dimensional  error  of  a  position,  as  applied  to 
the  special  case  of  range-a ziauth  fixes,  must  next  be  exam¬ 
ined.  Two  figures,  the  ellipse  and  the  circle,  are  used  to 
characterize  two-dimensional  error.  The  "error  diamond"  is 
also  sometimes  used  but  has  no  statistical  significance 
[Thomson,  1977].  The  error  ellipse  is  discussed  first  since 
it  is  the  aost  general  index  of  error.  Another  is  root  mean 
square  distance  (drml  ) ,  also  known  as  root  mean  square  error 
[Bowditch, 1977 ],  which  is  the  radius  of  a  circular  figure 
commonly  used  in  hydrography.  It  is  the  error  index  used 
for  NOS  positioning  standards  [Ombaoh,  1976],  A  second 
circular  figure,  known  as  circular  standard  error,  is  also 
examined  briefly  because  of  the  ease  of  converting  it  to 
circular  figures  which  have  different  probabilities. 

Por  clarification  of  the  issues  involved  in  this 
section,  the  following  assumptions  are  made. 

(i)  Only  random  errors  are  considered. 

(ii)  Errors  associated  with  each  LOP  are  normally 
distributed. 

<iii)  Errors  are  independent. 

(iv)  Errors  are  limited  to  the  two-dimensional  case. 

(v)  LOP's  are  straight  lines  at  their  point 
of  intersection* 
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These  are  the  saae  assumptions  made  by  Kaplan  (1980), 
except  that  Kaplan  allows  systematic  error  to  be  considered 
in  assumption  (i).  This  appears  to  be  an  oversight  since 
the  remainder  of  his  discussion,  from  which  this  section 
draws  heavily,  considers  only  random  error. 

1.  Concentric  and  Eccentric  Seoiatxy 

Before  proceeding  further  into  a  discussion  of  error 
figures,  it  is  necessary  to  examine  the  two  special  cases  of 
range-azimuth  positioning,  that  of  eccentric  and  concentric 
geometry.  Each  is  illustrated  in  Figures  2.2  and  2.3. 


In  actual  practice  the  geometry  used  is  often  eccen¬ 
tric,  but  the  concentric  assumption  is  made.  This  is 
because  it  is  usually  difficult  to  co-locate  both  theodolite 
and  ranging  equipment.  Hence,  they  are  offset  one  or  two 
meters  from  each  other.  It  should  be  noted  that  this 
assumption  will  introduce  a  systematic  error  in  all 
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Figure  2.3  Coscentric  Bange-lzimuth  Geometry. 

positions,  which  aust  be  eliminated  before  an  analysis  of 
random  errors  can  be  made.  The  systematic  error  is  often 
ignored  (by  use  of  the  concentric  assumption)  if  the  total 
uncompensated  error  is  within  the  tolerance  of  the  standards 
being  used.  An  algorithm  for  eliminating  this  error  is 
shown  in  Figure  2.4  and  given  by  eguation  2.6.  This  algo¬ 
rithm  assumes  that  the  size  of  d  in  Figure  2.4  is  small 
compared  to  r. 


Cad  COS  0 


(2.6) 


where: 

c  «  the  corrector  to  be  applied  to  r 
r  *  observed  range  to  the  vessel 

d  «  distance  between  theodolite  and  ranging  device 
f  >  the  angle  between  the  visual  LOP  and  the 
line  connecting  the  two  stations 
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Figure  2.4  Eccentric  Error  Compensation. 


2.  1M  atas.  Slii£§e 

Detailed  discussions  of  the  development  of  the  error 
ellipse  can  be  found  in  many  references,  especially  in 
Greenwalt  (1962)  and  in  Burt  (1966).  This  paper  will  only 
present  enough  background  to  apply  the  error  ellipse  concept 
to  two  LOP  range-azimuth  positioning.  The  error  ellipse 
formed  when  multiple  LOP  observations  are  made  is  not 
considered  here. 

h  range-azimuth  position  is  formed  by  the  intersec¬ 
tion  of  two  LDP's,  each  having  an  associated  standard  devia¬ 
tion.  By  applying  the  two-dimensional  normal  distribution 
to  the  errors,  elliptical  contours  of  equal  probability 
density  are  formed.  The  contours  center  on  the  intersection 
point  of  the  lines  of  position.  This  is  illustrated  in 
Figure  2.5,  and  shown  mathematically  by 
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where : 


v4  =  residual  (error)  in  the  direction  of 
vfc  *  residual  (error)  in  *  ha  direction  of  C £ 
*  length  of  the  semi-major  axis 
<3^*  length  of  the  seal-minor  axis 


It  can  be  shown  that 


V  X  K 

A  *  \  * 

-r  — r  =  K 

\  <rj 


(2.8) 


where 


—  2. 


[p  (v.,  Vb) 


<r*  <fb  2  rr 


(2.9) 


•'For  values  of  p  (v  ,  v  )  froa  0  to  oo  ,  a  faaily  of  equal 
probability  density  ellipses  are  foraai  with  axes  K  <&  and 
Kcrw"  [Greenwalt,  1962].  The  probability  density  function  in 
equation  2.7,  when  integrated  over  a  particular  area, 
becomes  the  probability  distribution  function.  This  yields 
the  probability  that  the  residuals  v*  and  vb  will  occur 
simultaneously  within  that  region.  This  probability  distri¬ 
bution  function  of  an  ellipse  is  given  by 


p  K  ,vb )  »  1  -  e 


(2.  10) 


The  solution  of  equation  2.8  for  different  values  of  K 
yields  different  probabilities.  For  exaaple,  for  39.35!? 
probability,  the  axes  of  the  ellipse  are  1.000  0^  and  1.000 
cfh  [Greenwalt,  1962  ].  In  other  words,  a  one-sigma  error 
ellipse  around  a  measured  position  indicates  a  39.35*  prob¬ 
ability  that  the  position  is  actually  within  that  ellipse. 

It  is  seen  in  Figure  2.5  that  the  standard  devia¬ 
tions^  and  Cl  of  the  measured  LOP's  are  not  the  same  as  the 
standard  deviations  Cl  and  <J£  of  the  error  ellipse.  A 
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Figure  2.5  Error  Ellipse  and  dr*t  . 


coordinate  transforeation  is  required  to  obtain  them.  This 
transformation  is  found  in  Heinzen  (1977)  and  will  not  be 
discussed  here.  Results  of  the  transformation  are  presented 
in  equations  2.11  and  2.12,  from  Bowditch  (1977). 


where 


0*  *  length  of  the  seei-major  axis 
*  length  of  the  seai-ainor  axis 
oi  »  standard  deviation  of  the  range  LOP 
Oi  =  standard  deviation  of  the  angle  LOP 
when  converted  to  distance  units 

1 9  ■  angle  of  intersection  of  LOP's 

For  range-azimuth  positioning,  C,  and  (Tx  are  not 
equal.  <T,  is  the  error  in  distance  aeasureaent,  and  it  is 
dependent  on  the  equipaent  used  for  ranging,  k  diagram  of  a 
range-azimuth  position  is  shown  in  Pigure  2.6. 


Flgare  2.6  &  Range- Axiaath  Position 


The  error  in  the  visual  LOP,  <TX  ,  is  a  function  of  the 
angular  error  of  the  theodolite  and  the  distance  to  the 
vessel,  where  r  is  the  distance  and  0^  is  the  angular  stan¬ 
dard  deviation  of  the  theodolite,  in  units  of  degrees.  This 
is  given  by  equation  2.13,  which  is  a  aodification  froa 
Heinzen  [1977].  Heinzen  uses  the  tarn  angular  resolution  in 
place  of  the  aore  correct 


*  S 7.2*0 


(2.  13) 


The  error  ellipse  concept  can  now  be  applied  to  the 
eccentric  geoaetry  by  using  equation  2.11,  2.12,  and  2.13 
directly.  In  the  concentric  case,  since  the  angle  of  inter 
section  , ,  is  always  ninety  degrees,  equations  2.11  and 
2.12  can  be  siaplified  to  equations  2.14  and  2.15. 


(2.  14) 


(2.  15) 


3.  S22S  sain  asms  Bis*  mss 

Boot  aean  square  distance  (drni{)  is  presented  here 
because  of  its  cosaon  use  in  hydrography.  It  is  not 
coaaonly  known  aaong  hydrographers  that  unlike  the  error 
ellipse,  drmthas  a  variable  probability  depending  on  the 
eccentricity  of  its  associated  error  ellipse,  and  ranges 
froa  68. 3*  to  63.2X,  as  shown  in  Figure  2.7  froa 
[Burt, 1977].  Eccentricity  is  defined  as  the  ratio  of  the 
seai-aincr  to  the  seai-aajor  axes  of  the  error  ellipse. 

Boot  aean  square  distance  is  also  called  Bean  Square 
Positional  Error  (BPSE)  by  Greenwalt  (1962),  who  recoaaends 
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1 


Figure  2.7  Variation  of  irms  Probability. 


that  this  index  not  be  used  because  of  its  variation  in 
probability. 

Soot  aean  square  distance  is  defined  in  Bovditch 
(1977)  as  equation  2.16.  An  alternative  fora  is  given  by 
Heinzen  (1977)  as  equation  2.17.  if  the  errors  are  assumed 
independent,  the  correlation  coefficient  is  zero  and 
equation  2.17  is  reduced  to  equation  2.18.  An  alternate 
method  of  arriving  at  equation  2.18  is  to  substitute  equa¬ 
tions  2.11  and  2.12  directly  into  equation  2.16. 


dr<**>  * 


(2.  18) 


If  concentric  geoeetry  is  assuaed,  the  angle /3  is  equal  to 
ninety  degrees  and  equation  2.18  is  reduced  to  equation 
2.19. 


(2.  19) 


Substituting  equation  2.13  into  equation  2.19, 
for*  for  range-aziauth  dKms  is  obtained. 


the  final 

(2.20) 


*•  Sir  gal ai  SlaMari  Susi 

Circular  standard  error  has  experienced  little  use 
in  hydrography  but  is  valuable  because  it  allows  easy 
conversion  between  circles  of  different  probability.  It  is 
derived  in  Graenwalt  (1962)  and  given  by  equation  2.21.  It 
should  be  noted  that  this  equation  is  only  an  approximation, 
although  a  very  good  one. 

Cl  -  0.5000  (  <T«  ♦  <Tb  )  (2.21) 

Circular  standard  error  has  a  probability  of  39.35*  for  a 
completely  circular  error  ellipse.  It  is  preferred  over 
drms  because  it  can  be  converted  to  other  circular  error 
indices  cf  different  probability  by  a  constant  conversion 
factor,  as  long  as  the  ratio  of  is  between  0.2  and 

1.0.  The  equation  for  circular  standard  error,  and  for 
other  circular  error  figures,  is  given  in  Table  I,  which  is 
taken  directly  froa  Greenwalt  (1962).  This  fable  gives  all 
error  indices  in  teras  of  either  <TC  or  the  error  ellipse 


TABLE  I 

Circular  Error  Foraulae 


Precision  Index 


Percentage 


Circular 

Standard  Error 

39- 35* 

Circular 

50% 

Probable  Error 

Circular  Map 

90% 

Accuracy  Standard 

Circular  Near- 

99-78% 

Certainty  Error 

(Three -five  sigma) 

Formula 


0.5000  (0x 

when  °min/° 


CMAS  =  1.0730  (ox 
when  u.Ja 


3  -5000  a 


seai-iajor  and  seai-iinor  axes  0^  and  <J£  .  Applying  the 
assuaptions  of  the  previous  section  on  drm{  to  this  case,  0~c 
for  the  concentric  range-aziauth  case  is  given  by  equation 
2.22. 


<n.  -  0.5000  57 


r 

7.296 


(2.22) 


C.  THE  EBBOB  OF  AN  INTERPOLATED  FIX 

The  error  in  an  interpolated  position  for  both  methods 
of  interpolation  used  in  range-azimuth  positioning  will  now 
be  derived.  Fundamental  to  the  derivation  is  an  under¬ 
standing  of  the  concept  of  error  propagation,  which  is 
explained  in  detail  in  Greenwalt  (1962)  and  Davis  (1981). 
Error  propagation  is  summarized  here  for  the  special  case  of 
range-azimuth  positioning.  It  must  be  noted  that  any  inter¬ 
polation  discussed  in  this  chapter  is  strictly  due  to  errors 
in  the  observed  positions  between  which  the  interpolation  is 
made.  Error  due  to  the  vessel  not  being  at  its  interpolated 
position  (due  to  steering  or  wind  and  sea  conditions)  will 
be  considered  in  later  chapters. 

i.  iB&asasla&laa  UastiUiis 

The  present  NOS  methods  of  interpolating  range- 
azimuth  fixes  are  of  two  types.  One  interpolates  both  the 
range  and  angle  between  two  observed  positions.  The  second 
is  used  when  actual  range  information  is  acquired  on  each 
sounding.  In  this  case  only  the  angle  is  interpolated,  and 
is  used  with  the  observed  range  to  compute  a  position.  In 
each  case  a  linear  interpolation  is  used  [Ehrhardt,  1979]. 
Algorithms  for  the  interpolated  value  of  the  range,  rx ,  and 
the  angle  6^,  are  given  in  equations  2.23  and  2.24. 

*  r.  *  CJ/R+1>  <rx-r,>  (2.23) 


-  e,  ♦  (j/K+i)  (e*-e, ) 


(2.24) 


where: 

Subscript  i  denotes  interpolated. 

Subscript  1  denotes  observed  position  number  1. 
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Subscript  2  denotes  observed  position  nuaber  2. 

K  denotes  the  nuaber  cf  interpolations, 
j  steps  froa  one  to  K. 

The  following  ezaaple  of  the  interpolation  process 
is  also  sketched  in  figure  2.8. 


Figure  2.8  Ezaaple  of  ingular  Interpolation. 


Let:  0,  *  251°  r,  *  1501m  j  *  1  to  2 

ez  *  244°  r*  »  1485m  K  *  2 

The  values  of  the  interpolated  angles  and  ranges  are 
coaputed  by  equations  2.24  and  2.23. 

8ai  »  251  ♦  (1/3)  (244-251)  -  251  -  2.3  *  248.7°  (2.25) 

0ai  ■  251  ♦  (2/3)  (244-251)  -  251  -  4.6  »  246.4° 


(2.  26) 


rv,  ■  1501  ♦  (1/3)  (1485-1501)  *  1501-5.3  *  1495.7m  (2.27) 


riz  «  1501  ♦  (2/3)  (1485-1501)  *  1501-10.7  *  1490 . 3m(2.  28) 


2.  JlISI  E£OB&aa&g& 


"Error  propagation  is  better  termed  the  propagation 
of  variances  and  covariances"  [Davis,  1981].  The  following 
paragraphs  are  a  general  derivation  of  error  propagation, 
and  will  be  applied  below.  In  reading  this  section,  the 
terns  x  and  y  are  general,  yet  tha  reader  should  remember 
that  they  will  be  applied  specifically  to  range-azimuth 
interpolation.  Let  y  be  a  set  of  quantities  each  of  which 
is  a  function  of  another  set  of  random  variables  x.  The 
random  variables  in  our  application  are  the  observed  x 
given  in  equation  2.  1,  and  y  is  the  interpolated  or  ru 
given  in  equations  2.24  or  2.23. 

The  covariance  matrix  yy  is  given  by  the  matrix 
equation  2.29,  where  Jy„  is  called  the  Jacobian  matrix  and 
is  given  in  equation  2.30. 


E 


yy 


(2.29) 


(2.30) 
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is  the  transpose  of  Jy*  ,  and  2jXX  is  another  covariance 
matrix  given  by  : 


(2.31) 


The  covariance  matrix  in  equation  2.31  is  written 
for  the  general  case  of  correlated  random  variables.  If  y 
is  a  single  quantity  rather  than  a  set  of  quantities,  and  if 
the  variables  are  assumed  uncorrelated,  equation  2.29 
reduces  to: 


(2.32) 


By  carrying  out  the  multiplication  in  equation 
expanded  form  is  given  in  equation  2.33. 


2.32,  the 


(2.  33) 


If  y  is  a  linear  function  of  random  variables,  the 
partial  derivative  terms  in  equation  2.33  become  constants. 
Thus  the  matrix  equation  2.29  becomes  equation  2.34,  where 
a,b,c,  and  d  are  constants. 


*  i 

a  <r, 


I  *  * 

b  <rx 


+  Cf5  +  —  - 


d  <r* 


(2.34) 
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Mow  1st  us  apply  equation  2.34  to  the  range-azimu-h 
case.  The  assuaptions  made  above  are: 

(i)  y  is  a  single  quantity 

(ii)  The  randoa  variables  x  are  assumed 
uncorrelated. 

(iii)  y  is  a  linear  function  of  x 

If  these  are  applied  to  the  two  intarpolation  algorithms, 
equations  2.24  and  2.23,  it  can  be  seen  that: 

(i)  0;  and  are  single  quantities. 

(ii)  The  randoa  variables  8,  ,  d2  and 
r,  ,  r4  are  assuaed  uncorrelated. 

(iii)  The  two  interpolation  algorithms  are  linear. 


From  the  foregoing  general  discussion  the  interpo¬ 
lated  variances  of  equations  2. 23  and  2.24  can  be  given  as: 


*•♦(£-)(**- ,r") 

but, 

*  <J> 


(2.35) 


(2.36) 


(2.  37) 
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B 


and: 


<£,  *  <£*  *  <*0  (2»  38) 

By  substituting  equations  2.37  and  2.33  into  equations  2.35 
and  2.36,  we  have: 


a 


cs*  = 


and 


del 


(2.39) 


(2. 40) 


This  error  propagation  applies  only  to  vessels 
moving  in  an  arc.  Por  a  vessel  moving  in  a  straight  line 
the  angular  interpolation  algorithm  is  not  linear  and  the 
partial  derivative  terms  analagous  to  those  in  aquation  2.33 
are  not  constants.  The  error  ellipse  for  an  interpolated 
position  can  now  be  formed.  The  ellipse  is  seen  to  be  the 
same  as  for  an  observed  position  (figure  2.5)  since  equa¬ 
tions  2.39  and  2.40  show  that  errors  in  interpolated 
distance  and  angle  are  the  same  as  for  an  observed  distance 
and  angle. 
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III.  limimi  5I5IS*  m  IKLSIIKAIIQI 


The  experimental  work  performed  in  this  investigation 
can  be  conceptually  divided  into  two  parts,  although  both 
parts  were  accomplished  simultaneously.  Part  one  involved 
pointing  error  and  variance  determination  for  the  two  theo¬ 
dolites,  and  part  two  investigated  the  accuracy  of 
interpolated  positions. 

A.  FIELD  WORK 

The  actual  field  work  was  representative  of  a  typical 
range**aziauth  survey.  Two  full  days  (8  and  15  April,  1983) 
were  required  to  obtain  443  position  fixes  and  over  2500 
Aztrac  angles.  The  experiment  took  place  in  southern 
Monterey  Bay  near  the  Monterey  Harbor  Coast  Guard  Pier,  as 
shown  in  Figure  3.1.  The  vessel  used  was  a  chartered 
36-foot  Uniflize  with  a  fiberglass  hull  and  twin  engines, 
and  its  operator  had  about  two  months'  hydrographic  survey 
experience,  including  steering  range  arcs. 

For  each  position  of  the  vessel,  six  lines  of  position 
were  observed.  Three  Wild  T-2  theodolites  located  at 
stations  HOSSEL,  SOFAR,  and  USE  MON  were  used  to  obtain  the 
"reference"  or  best  estimate  positions  of  the  vessel.  Two 
additional  test  theodolites,  a  Wild  T-2  and  an  Odom  Aztrac, 
made  observations  free  stations  T2  and  AZTRAC.  Finally,  a 
Del  Norte  Trisponder  (model  R04)  provided  a  distance  LOP  to 
the  vessel  from  station  GEOCEIVER .  The  reference  positions 
were  obtained  at  one  minute  intervals,  and  five  Aztrac 
angles  were  observed  between  each  of  these  positions  for 
later  use  in  evaluating  interpolation  methods.  Figure  3.2 
shows  the  sources  of  the  various  LOP's  to  the  vessel.  Zt 
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Figure  3.1  Sketch  of  the  Survey  tree. 


represents  an  enlargeaent  of  the  dashed  circle  area  shown  in 
Figure  3.1. 

Geodetic  control  for  the  reference  positions  consisted 
entirely  of  aonuaented  third-order  stations.  Station 
GEOCEIVER,  although  not  a  published  station,  was  located  to 
third  order  specifications  by  Hr.  Williaa  Anderson  of 
NAVOCEANO  in  1982,  and  the  other  two  stations  were  located 
as  eccentrics  of  this  station.  Each  was  less  than  one  seter 
fros  GEOCEIVER.  GEOCEIVER  is  a  aonueented  station,  while 
Aztrac  and  T2  are  aarked  by  aasonry  nails  driven  into  the 
concrete  pier.  Positions  for  Aztrac  and  T2  were  coaputed  on 
an  HP-9815  coaputer  using  the  NGS  geodetic  direct  prograa. 
The  initial  pointing  by  both  instruaents  was  to  station  OSE 


HON.  Initial  pointings  at  all  stations  were  made  to  third- 
order  stations  at  least  500  meters  distant.  Geodetic  posi¬ 
tions  of  all  control  stations  is  given  in  Appendix  B. 

Underway  operations  were  very  similar  to  a  nonautomated 
range-azimuth  survey,  except  that  five  theodolites,  rather 
than  one,  were  trained  on  the  boat  for  each  fix.  The  survey 
boat  steered  along  the  appropriate  range  arc,  and  fix  marks 
were  given  over  voice  radio  to  the  observers  on  shore.  All 
T-2  angles  were  obtained  once  each  minute,  and  Aztrac  angles 
were  recorded  every  ten  seconds.  The  sequence  of  events  for 
each  reference  position  is  given  in  Table  II. 

This  process  continued  throughout  the  two-day  field 
operation.  Breaks  in  data  collection  occurred  at  the  ends 
of  each  range  arc,  and  also  when  theodolite  observers  were 
rotated.  The  monotony  of  the  events  shown  in  Table  II  was 
sufficent  to  approximate  the  monotony  (with  its  associated 
observer  difficulties)  of  an  actual  survey. 
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TABLE  II 

Data  Acquisition  Sequence  of  Events 


time  prior 
to  fix 


20  sec 
15  sec 
10  sec 
5  sec 
0  sec 


actio  n 

depress  Aztrac  manual  fix  button 

"15  second  standby"  over  voice  radio 

depress  Aztrac  manual  fix  button 

"standby"  over  voice  radio 

"■art"  over  voice  radio, 
depress  Aztrac  manual  fix  button, 
T-z  operators  mate  observations 


Observers  at  stations  MUSSEL  and  USE  MON  were  required 
to  record  angles  as  well  as  to  observe.  T-2  observers  at 
stations  T2  and  SOFAS  were  provided  with  separate  recorders, 
because  at  times  the  high  angular  speeds  of  the  vessel  at 
these  stations  (due  to  closeness  of  the  boat)  made  observing 
and  recording  difficult  for  one  person. 

The  Trisponder  was  calibrated  using  the  standard  NOS 
method  over  a  geodetic  baseline  (GEOCEIVER  to  USE  HON)  and 
no  systematic  errors  were  observed.  The  geodetic  baseline 
was  deterained  by  coaputing  an  inverse  distance  between 
stations  USE  MON  and  GEOCEIVER  using  the  HP9815  computer  and 
NGS  geodetic  software.  The  Trisponder  was  reported  by  Odom 
Offshore  Surveys  to  have  a  standard  deviation  of  its  ranging 
error  of  1.0  aeters.  The  master  unit  and  antenna  were 
mounted  at  the  highest  part  of  the  boat,  about  one  meter 
above  and  aft  of  a  radar  antenna  enclosure,  which  was  about 
one  meter  in  diameter  and  0.4  meter  high.  The  master  unit 
was  a  cube  about  0.3  meter  on  a  side  and  was  covered  with 


Tr "  *  'xem 


‘vswmww.i  v  m  jjpjjjjjippgj 


green  signal  cloth.  This  unit  vas  the  target  to  which 
angles  were  observed.  The  Trispondar  Distance  Measuring 
Unit  (DMU)  was  nounted  next  to  the  steering  station  inside 
the  boat. 

Data  logging  for  both  Trisponder  and  Aztrac  was  accom¬ 
plished  by  an  Odom  Nartrace  computer  system.  Both  the 
Trisponder  DMU  and  Aztrac  receiver  were  interfaced  to  this 
unit.  The  Navtrace  computer  is  programmed  to  automatically 
log  fix  data  not  on  intervals  of  equal  time,  but  on  equal 
distance  intervals  from  a  reference  line.  Por  this  reason, 
all  fixes  were  logged  by  using  the  manual  fix  feature  of  the 
computer.  This  simply  caused  the  computer  to  log  distance 
and  angle  each  time  a  button  was  manually  pressed.  Timing 
for  fixes  was  provided  by  an  HOS  standard  sounding  clock. 
This  is  a  mechanical  clock  with  a  buzzer  set  to  ring  every 
ten  seconds. 

The  weather  during  both  days  of  field  work  was  good  zo 
fair  for  survey  work.  Visibility  was  good  at  all  times,  and 
winds  were  calm  on  each  morning  of  operations.  Afternoon 
northwest  winds  were  a  maximum  of  about  15  knots  on  both 
days,  which  produced  two  to  three  foot  seas  in  the  offshore 
part  cf  the  operating  area. 

B.  THEODOLITE  POINTING  EBROB 

The  pointing  error,  xi#  is  defined  as  the  difference 
between  an  observed  and  computed  angle  given  by 

\  M  %  -  9t  (3.1) 

and  the  mean  pointing  error,  7,  is  therefore 


A*t 
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where  n  is  the  number  of  observations,  &a  is  the  observed 
azimuth  and  Bc  is  the  computed  azimuth  from  the  theodolite 
position  to  the  reference  position  of  the  vessel.  The  stan¬ 
dard  deviation,  <f9  ,  of  a  set  of  pointing  errors  can  be 
expressed  maths mat ic ally  as 

<5  -  <3-3> 

V  A*' 

This  is  the  usual  definition  in  most  statistics  texts  and 
agrees  with  the  discussion  in  section  II. A.  of  this  paper. 

The  ccaputed  azimuth,  9C  ,  was  determined  by  computing  a 
geodetic  forward  azimuth  from  either  test  theodolite  posi¬ 
tion  to  the  reference  position  of  the  vessel.  There  is  some 
potential  error  in  0C  due  to  an  uncertainty  in  the  reference 
position  of  the  vessel.  The  reference  position  of  the 
vessel  was  determined  by  a  least  squares  adjustment  of  the 
three  LOP's  from  the  three  Wild  T-2  theodolites.  A 
by-product  of  this  adjustment  is  an  error  ellipse  for  each 
position.  The  angle  subtending  this  error  ellipse  from  each 
test  theodolite  is  the  error  in  9e  .  This  is  illustrated  in 
Pigure  3.3. 


Figure  3.3  Uncertainty  of  an  Observed  Error 


The  least  squares  eethod  was  developed  in  1794  by  Gauss. 
It  is  based  on  the  ccncept  that  as  soon  as  redundant  obser¬ 
vations  are  present  there  appear  discrepancies  among  the 
observations.  In  order  to  eliminate  these  discrepancies  a 
residual,  v ,  must  be  added  to  each  observation.  In  the 
method  of  least  squares  the  residuals  are  determined  such 
that  the  sum  of  their  squares  becomes  a  minimum,  provided 
all  observations  have  the  same  accuracy  [Mueller,  1979].  It 
is  assumed  that  the  three  LOP»s  used  in  this  adjustment  do 
in  fact  have  the  same  accuracy.  Further  information  on  the 
method  of  least  squares  as  applied  to  hydrography  is  given 
by  Kaplan  (1980) . 

Reference  positions  were  calculated  by  the  least-squares 
FORTRAN  program  &ZLSQ2,  written  by  the  author,  which  is  a 
modification  of  program  SILVA1  [Silva,  1979].  The  program 
generated  plane  and  geodetic  coordinates  for  each  position 
fix,  as  well  as  the  lengths  of  the  semi-major  and  semi- minor 
axes  of  its  error  ellipse.  The  angle  mada  by  the  semi-major 
axis  and  the  x-axis  was  also  computed. 

C.  INTERPOLATION  ALGORITHM  EVAL0ATION 

The  present  NOS  method  of  observing  azimuths  to  the 
sounding  vessel  is  totally  nonautomared,  as  explained  in 
section  A.  This  precludes  recording  an  azimuth  for  each 
sounding  because  of  the  speed  with  which  soundings  are 
taken.  Consequently  interpolation  algorithms  are  employed 
to  plct  soundings  between  observed  positions.  Two  algo¬ 
rithms  are  used:  one  interpolates  only  the  azimuth  and  uses 
an  observed  distance  to  the  vessel,  while  the  other  interpo¬ 
lates  both  azimuth  and  distance.  Figures  3.4  and  3.5 
illustrate  the  situation. 
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Figure  3.4  Interpolation  of  Angle  Only. 
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Figure  3.5  Interpolation  of  Angle  and  Distance. 


The  nethod  used  to  evaluate  these  interpolation  aethods 
is  sinple.  Actual  distances  and  angles  to  the  vessel  were 


recorded  every  ten  seconds,  using  the  telemetering  theo¬ 
dolite  and  associated  computer  logging  system,  A  ncr.auto- 
mated  system  vas  then  simulated  by  interpolating  between 
positions  obtained  every  minute.  Thus  two  data  sets  were 
produced,  an  interpolated  set  and  a  corresponding  set  of 
actually  observed  positions.  These  two  sets  of  positions 
were  then  compared  and  their  differences  examined.  Both  NOS 
interpolation  algorithms  were  evaluated  in  this  manner.  An 
analysis  cf  the  experimental  results  is  given  in  Chapter  IV. 

D.  CHOICE  OF  EXPBHIHBNT1L  CONDITIONS 

Any  experiment  must  be  performed  under  conditions 
similar  to  those  under  which  the  results  will  be  applied. 
This  section  will  describe  the  specific  error  sources  to  be 
examined  under  range-azimuth  conditions.  By  carefully 
choosing  the  experimental  conditions,  these  error  sources 
can  be  brought  into  the  foreground  for  examination,  while 
all  ether  sources  of  error  can  be  kept  in  the  background. 

It  is  necessary  to  introduce  some  standard  statistical 
terminology  to  help  understand  the  experiment  design.  The 
classical  experimental  method  studies  the  effect  of  only  one 
variable.  That  is,  it  holds  all  effects  but  one  to  be 
constant,  and  varies  that  single  one  systematically,  in  the 
case  of  this  thesis,  one  factor  (a  theodolite)  was  varied 
systematically  by  introducing  two  levels  (Aztrac  and  T-2)  of 
the  factor.  It  is  desired  to  find  the  effect  of  that  factor 
(the  theodolite)  on  the  error  of  a  rang  8-azimuth  position. 

It  is  known,  however,  that  the  observer  also  has  a  consider¬ 
able  effect  on  the  error  of  a  position.  Thus  another  factor 
(the  observer)  is  introduced  into  the  experiaent.  Four 
levels  (four  different  individual  observers)  of  this  factor 
were  selected. 
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The  following  four  observers,  who  are  attached  to  the 
Naval  Postgraduate  School,  constitute  the  four  levels  of  the 
observer  factor. 

LCDS  Gerald  B.  Hills,  NO  A  A.  Instructor  in 
Hydrography.  Ten  years  experience  in 
geodetic  and  hydrographic  surveying. 

LT  Maureen  Kenny,  NOAA.  Student  in 
hydrography.  2  years  field  experience 
aboard  the  NOAA  Ship  Davidson,  primarily 
in  Alaskan  and  Hast  Coast  waters. 

LT  Hary  C.  Schoaaker,  NOAA.  Student  in 
hydrography.  4  years  field  experience 
aboard  NOAA  Ship  Davidson,  and  on  NGS 
horizontal  control  and  leveling  parties. 

Hr.  James  R.  Cherry,  Supervisory 
Geodesist,  Naval  Postgraduate  School.  23 
years  field  experience  in  hydrography  and 
geodesy,  with  the  Naval  Oceanographic 
Office. 

Although  additional  observers  would  have  been  desired,  none 
were  available  who  had  any  experience  with  theodolites. 

Aside  from  the  factors  to  be  investigated,  it  was  recog¬ 
nized  that  there  were  background  conditions  that  affected 
the  results  of  the  experiment  (Crow,  1955].  Some  of  these 
were  taken  into  account  explicitly  in  the  design.  The 
influence  of  all  others  was  minimized  by  scheduling  the 
experiment  such  that  each  combination  of  instrument  and 
observer  was  evaluated  in  random  order,  thereby  randomizing 
the  effects  of  these  other  conditions. 


Blunders  and  sone  systematic  errors  were  explicitly 
accounted  for  in  the  design.  Blunders  were  identified  visu¬ 
ally  by  their  large  magnitude  and  were  simply  deleted  from 
the  data  set  using  the  editor  on  the  UPS  computer.  The 
& ztrac  instrument  essentially  eliminates  all  blunders  due  to 
its  automatic  data  logging  feature. 

Systematic  errors  integral  to  the  theodolites,  such  as 
collimation  and  eccentricity,  were  not  accounted  for,  since 
they  were  very  small  compared  to  the  size  of  the  errors 
under  investigation.  One  systematic  error  that  could  be 
quite  large  is  the  initial  pointing  error,  which  is  due 
largely  to  the  observer.  The  NOS  Hydrographic  Manual 
[Ombach,  1976]  requires  the  initial  pointing  to  be  "accurate 
to  within  ±30  seconds  of  arc".  That  is,  the  difference 
between  beginning  and  ending  pointings  is  not  to  exceed  one 
minute.  For  this  experiment,  the  mean  of  the  beginning  and 
ending  pointings  to  the  initial  azimuth  was  algebraically 
subtracted  from  all  observed  angles  for  that  set.  It  could 
be  argued  that  a  more  accurate  way  of  removing  this  system¬ 
atic  error  would  be  to  prorate  the  difference  in  initial 
readings  between  the  beginning  and  ending  pointings. 

However,  the  former  method  is  often  used  by  NOS  hydrographic 
survey  units. 

One  background  condition,  the  angular  speed  of  the 
vessel,  was  explicitly  taken  into  account  in  the  design  by 
dividing  the  experiment  into  groups.  Six  separate  subsets 
of  the  experiment  were  created  by  measuring  pointing  error 
at  six  different  angular  speeds.  These  speeds  were  chosen 
to  be  closely  representative  of  speeds  found  under  actual 
surveying  conditions.  To  this  end,  integral  values  of 
angular  speed  were  not  used,  but  rather  integral  distances 
from  the  ve^al  to  the  theodolite  station.  This  corresponds 
to  actual  practice  in  most  nonautomated  or  semiautomated 
surveys.  The  vessel  maintains  constant  engine  speed,  and  is 
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navigated  in  a  circle  such  that  a  constant  distance  value  is 
displayed  in  the  ranging  eguipaent.  Par  this  experiment, 
distances  of  300,  500,  700,  1000,  1500,  3000  meters  were 
used.  Using  the  approximate  vessel  speed  of  2  aeters  per 
second  (4  knots)  ,  the  angular  speeds  in  minutes  per  second 
corresponding  to  these  distances  were  22.9,  13.8,  9.8,  6.9, 
4.6,  and  2.3,  respectively. 

Another  background  condition  that  was  explicitly 
accounted  for  in  the  design  was  the  axperience  of  the 
observers.  For  a  completely  general  investigation,  one 
would  desire  to  evaluate  the  error  using  observers  of  widely 
varying  experience.  Only  experienced  observers  were  used  in 
order  to  conserve  resources  of  time  and  money,  although  none 
of  them  were  experienced  in  using  the  new  Aztrac  instrument. 
All  observers  except  LCDS  Bills  had  acted  as  observer  for 
range-azimuth  hydrography  within  tha  past  year  using  the 
Bild  T-2. 

There  were  other  background  conditions  that  affected  the 
results  of  the  experiment,  including  weather  and  lighting 
conditions,  and  observer  fatigue.  Other  subtle  factors  may 
have  also  contributed  to  the  error.  In  nonautoaated 
systems,  there  might  have  been  a  time  lag  by  the  radio  oper¬ 
ator  in  the  vessel  as  the  fix  mark  was  relayed  to  the  theo¬ 
dolite  operator  ashore.  The  design  of  the  experiment  could 
not  explicitly  account  for  all  these  conditions,  so  several 
steps  were  taken  to  randomize  the  order  of  observations  for 
each  combination  of  instrument  and  observer. 

It  was  assumed  that  the  theodolites  used  were  randomly 
drawn  from  the  entire  population  of  T-2  and  Aztrac  instru¬ 
ments.  This  is  a  fairly  good  assumption  for  the  T-2,  since 
the  particular  unit  tested  is  one  of  several  maintained  by 
the  Naval  Postgraduate  School,  and  has  seen  several  years  of 
service.  The  assumption  for  the  Aztrac  theodolite  is  not  as 
good,  since  it  is  one  of  only  five  in  existence  at  this 
time.  It  was  provided  by  the  manufacturer  for  testing. 
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Theodolite  observers  were  required  to  be  paired,  sines 
the  experimental  theodolites  observed  simultaneous  angles  to 
the  vessel,  A  random  selection  of  observer  pairs  was  accom¬ 
plished  by  consulting  a  random  number  table  in  Wonnacott 
(1977).  The  order  of  observation  for  the  six  different 
range  arcs  was  important  and  required  randomization. 

Failure  to  do  so  could  allow  the  error  associated  with  each 
to  be  influenced  by  time  varying  conditions.  Examples  of 
these  are  the  changing  effect  of  sun  glare  on  the  instru¬ 
ment,  observer  fatigue,  the  effect  of  repetition  acting  to 
decrease  error,  or  increasing  afternoon  winds  disturbing  the 
theodolite.  Randomization  was  again  accomplished  with  the 
aid  of  a  random  number  table.  A  random  number  was  assigned 
to  each  range  arc,  and  the  order  of  observation  was  estab¬ 
lished  by  selecting  the  arcs  from  the  lowest  to  highest 
number. 
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A.  DATA  PROCESSING  SYSTEM 

All  data  acquired  for  this  experiment  wer  3  manually 
entered  into  the  NPS  coaputer  facility.  A  data  processing 
systea  was  designed  and  implemented  to  log  the  raw  data  and 
perform  the  necessary  computations  to  arrive  at  the  finished 
form  given  in  this  chapter.  This  system  was  divided  into 
numerous  subsystems  to  entar  data,  determine  the  pointing 
error,  compute  interpolated  positions,  compute  means  and 
variances,  test  for  randomness,  and  compute  analysis  of 
variance  (ANOVA)  statistics.  ANOVA  is  a  standard  statis¬ 
tical  technique  used  in  testing  for  a  difference  between  the 
effects  of  two  or  more  factors. 

The  data  processing  system  was  designed  and  programmed 
by  the  author.  As  mentioned  previously,  the  NGS  geodetic 
inverse  and  direct  subroutines  (Pfeifer,  1975  ],  as  well  as 
the  least-squares  adjustment  program  [Silva,  1979  ],  were 
adapted  for  use  here.  The  ANOVA  computations  were  performed 
by  the  Statistical  Package  for  Social  Sciences  (SPSS)  on  the 
Naval  Postgraduate  School  mainframe  computer  (Hull,  198  1]. 
The  pointing  error  and  interpolation  subsystems  were 
designed  as  described  in  Chapter  III. 

B.  POINTING  BBROR  DETERMINATION 

The  pointing  error  standard  deviation  for  each  theo¬ 
dolite  needs  to  be  quantified,  and  a  determination  must  be 
made  as  to  whether  there  is  a  statistically  significant 
difference  between  the  two  instruments.  Pointing  error  for 
each  combination  of  instrument  and  observer  was  determined 
by  the  methods  discussed  in  Chapter  III.  Standard  deviation 
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of  the  pointing  error  was  computed  using  equation  3.3,  for 
each  combination  of  instrument  and  observer,  at  each  range 
arc,  and  are  given  in  Appendix  A.  This  entire  procedure  was 
performed  first  for  the  T-2  angles  observed  to  the  nearest 
second  of  arc,  which  is  the  maximum  resolution  for  this 
instrument,  and  second,  for  T-2  angles  rounded  to  the 
nearest  minute  of  arc,  which  is  specifically  allowed  by  the 
NOS  Hydrographic  Manual  [Umbach,  1976  ], 

Since  eight  different  instrument-observer  combinations 
exist,  the  sample  means  and  standard  deviations  for  each 
combination  were  slightly  different  because  each  sample  is 
only  an  estimate  of  the  mean,  ,  and  standard  deviation,  , 
for  the  entire  population.  The  pooled  standard  deviation  is 
the  best  estimate  of  for  this  case  of  multiple  samples 
because  it  takes  the  differences  among  sample  means  into 
account.  If  an  overall  standard  deviation  is  computed  using 
equation  3.3,  the  population  standard  deviation  will  be 
overestimated  because  of  the  differences  among  sample  means. 
The  pooled  standard  deviation  is  mathematically  expressed  by 
equation  4.1  ^crow,  1955  ],  [Box,  1978  ], 


>  (*>-Qsc  + . +(n*-i)*K 


n, 


- +  K 


where: 

k  *  total  number  of  observer-instrument 
combi  na  tions 

n  *  number  of  observations  for  each  observer- 
instrument  combination 
s  *  sample  variance  for  each  observer- 
instrument  combination 


Thus  a  pooled  standard  deviation  was  computed  for  the  eight 
samples  available  at  each  different  range  arc.  These  are 
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given  in  Table  III  for  the  Aztrac  and  both  the  rounded  and 
unrounded  T-2  data.  This  table  gives  angular  speeds  for  the 
vessel  in  units  of  minutes  of  arc  per  seconds  of  tiae.  Note 
that  the  pointing  error  of  the  Aztrac  is  twice  as  large  as 
that  of  the  T-2  and  that  the  rounded  T-2  pointing  error  is 
only  slightly  greater  than  the  unrounded. 


TABLE  III 

Pointing  Error  Standard  Deviation  (pooled  estiaates) 


Range 

Arc 

(a) 

Approx. 

Angular 

Speed 

(ain/sec) 

Aztrac 

a  (sec) 

T-2 

(unrounded) 

i  (sec) 

T-2 

(rounded) 

a  (sec) 

i  <71 

_ Oe _ 

<31 

<r. 

300 

22.9 

I  3.33 

(2  290) 

1.  26 

(868) 

1.31 

(902) 

500 

13.8 

|  3.56 

(1470) 

1.72 

(712) 

1.73 

(715) 

700 

9.8 

I  2.10 

| 

(6  19) 

1.30 

(382) 

1.31 

(3  86) 

1000 

6.9 

i  2.97 

(6  13) 

1.08 

(225) 

1.09 

(226) 

1500 

4.6 

1  4.25 

(584) 

0.92 

(126) 

0.93 

(128) 

3000 

2.  3 

1  2.06 

(142) 

1.41 

(  97) 

1.43 

(  99) 

The  technique  of  analysis  of  variance  (ANOV A)  is  used  tc 
deteraine  whether  a  statistically  significant  difference 
exists  between  the  two  theodolites.  The  ANOV A  technique 
allows  the  partitioning  of  overall  pointing  error  variance 
into  portions  caused  by  each  factor  (observer  and  theo¬ 
dolite)  ,  by  interaction,  and  by  experiaental  error. 
Interaction  exists  if  the  variance  for  a  particular  combina¬ 
tion  of  instruaent  and  observer  is  greater  than  the  variance 
for  any  other  such  coabination.  The  experiaental  error  is 


primarily  due  to  the  conditions  described  in  Chapter  III, 
concerning  the  reference  positions  for  the  vessel.  This 
error  is  a  neasure  of  the  precision  of  the  experiment,  where 
precision  is  defined  as  the  closeness  with  which  repeated 
measurements  made  under  similar  conditions  are  grouped 
together  [Green wait,  1971].  Experimental  error  is  assumed 
to  have  a  population  mean  of  zero.  Further  details  of  the 
ANOV k  method  can  be  found  in  many  statistics  textbooks,  such 
as  Wonnacott  (1977),  Box  (1  978),  Crow  (1955),  and  Halpole 
(1978).  This  discussion  was  taken  primarily  from  Crow 
(1955). 

The  following  assumptions  must  be  made  when  using  the 
analysis  of  variance  technique. 

(i)  Observations  are  random. 

(ii)  Means  and  variances  are  additive,  as  given 
in  the  mathematical  model  below. 

(iii)  Experimental  errors  are  independent. 

(iv)  Variances  of  the  experimental  errors 
are  equal. 

(v)  Distribution  of  the  experimental  errors 
is  normal. 

k  mathematical  model  can  cow  be  given,  using  these 
assumptions,  which  specifies  the  total  effects  of  the  vari¬ 
ances  acting  on  a  particular  observation,  x^+. 

V**  a,  ♦  bj  .  <«.2) 


where: 

■  overall  mean  for  all  observations 
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a^,  *  effect  of  the  instrument  factor  at  leva!  i 
b^  *  effect  of  the  observer  factor  at  level  j 

»  effect  of  the  interaction  of  instrument  and 
observer  at  level  i  and  j,  respectively . 

e^t  »  random  effect  caused  by  the  variance  of  the 

J 

experimental  error, 
i  a  1,2 
j  *  1,2, 3, 4 
t  3  1,2,3,  ...  ,  n 

n  «  the  number  of  observations  for  a 
particular  observer-instrument 
combination  (usually  15) 

It  should  be  emphasized  that  the  a,  b,  i  ,  and  e  of  equa¬ 
tion  4.2  are  not  actual  variances,  but  a  realization  of  the 
effect  of  those  variances  on  a  particular  observation  . 
The  variances  associated  with  observer,  instrument,  and 
interaction  are  computed  by  the  AN07A  procedure,  and  these 
form  the  basis  of  the  test  for  differences.  This  test  is 
called  the  P-test  in  honor  of  Sir  Ronald  A.  Pisher 
[Wonnacott,  1977  ],  and  is  based  on  a  ratio  of  variances.  To 
test  for  a  difference  between  instruments,  we  form  a  ratio 
with  the  variance  among  instruments  in  the  numerator,  and  a 
denominator  composed  of  an  estimate  of  the  variance  of  the 
experimental  error.  In  terms  of  our  mathematical  model,  a 
(the  numerator)  is  being  compared  with  eAit  (the  denomi- 
nator).  norm  simply  stated,  it  is  a  comparison  of  the 
precision  of  the  instruaent  (the  numerator),  with  the  preci¬ 
sion  of  the  experiment  (the  denominator) .  This  ratio,  F,  is 
then  compared  to  a  ratio  ^  , which  is  computed  for  for  a 
particular  confidence  level  from  the  F-distribution  function 
given  as  equation  4.3  [Crow,  1955  ]. 
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where:  F  *  the  ratio  discussed  above 

f,  *  the  number  of  degrees  of  freedom 
in  the  numerator  of  F 

fx  =  the  number  of  degrees  of  freedom 
in  the  denominator  of  P 

A  precise  hypothesis  must  now  ba  stated  that  can  be 
tested  by  the  P-test.  'Walpole  (1974)  states  that  "a  statis¬ 
tical  hypothesis  is  an  assumption  or  statement,  which  may  or 
may  not  be  true,  concerning  one  or  more  populations." 
Experiments  are  designed  to  test  hypotheses,  and  "the  rejec¬ 
tion  of  an  hypothesis  is  to  conclude  that  it  is  false,  while 
acceptance  merely  implies  that  we  have  no  reason  to  believe 
otherwise"  [Walpole,  1978].  A  null  hypothesis  is  an 
"initial  hypothesis,  or  one  we  hops  to  reject"  [Crow,  1955], 
and  is  usually  stated  in  terms  of  an  assumption  of  no 
difference  between  the  effects  to  be  investigated  by  the 
experiment.  This  experiment  uses  three  null  hypotheses, 
which  are: 

(i)  There  is  no  difference  between  observers, 
that  is,  the  ratio  F,  for  observers,  is 
small  compared  to  F* . 

(ii)  There  is  no  difference  between  instruments, 
that  is,  the  ratio  F,  for  instruments,  is 
small  compared  to  F* . 

(iii)  There  is  no  interaction,  that  is,  the  ratio 
F  for  interaction,  is  small  compared  to  P*. 
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These  three  null  hypotheses  can  be  accepted  or  rejected 
on  the  basis  of  the  P-test  and  the  ANOVA  procedure.  If  any 
are  to  be  rejected,  the  value  of  the  ratio  P  eust  exceed 
that  of  the  critical  P*  for  that  confidence  level  (95%)  .  .\ 

suaaary  of  A  NOV  A  results  is  given  in  Table  IV.  For  each 
range  arc,  the  coaputed  values  of  F  and  the  corresponding 
critical  value  needed  to  reject  the  null  hypothesis  are 
given.  The  rightaost  coluan  of  the  areas  labelled  as 
rounded  or  unrounded  data  indicate  acceptance  or  rejection 
of  the  null  hypothesis.  It  can  be  seen  froa  this  table  that 
in  no  case  could  the  null  hypothesis  be  rejected  for  either 
instruaent  or  observer.  In  other  words,  it  may  be  said  with 
95*  confidence  that  there  is  no  reason  to  believe  there  is  a 
difference  between  the  four  observers  or  between  the  two 
instruaents.  It  should  be  noted  that  data  for  the  1500 
aeter  arc  indicate  a  rejection  of  the  null  hypothesis  for 
interaction.  According  to  Crow  (1955),  a  significant  inter¬ 
action  usually  occurs  because  unrandoaized  background  condi¬ 
tions  are  present.  There  is  littla  apparent  reason  why  this 
interaction  should  occur,  other  than  that  there  were  some 
unidentified,  tiae  varying  conditions  affecting  the  measure- 
aents.  The  1500  and  the  1000  aeter  range  arc  were  only 
observed  on  April  8,  while  all  the  other  range  arcs  were 
investigated  on  both  days  of  the  experiment. 

A  final  assumption  of  the  ANOVA  technique  is  that  the 
observations  within  the  eight  combinations  of  instrument  and 
observer  are  randomly  drawn  froa  their  populations.  Several 
tests  are  available  to  deteraine  if  a  particular  sample  is 
randoa.  The  test  used  here  was  the  Run  Test,  as  given  in 
Crow  (1955)  ,  and  a  confidence  level  of  95*  was  used.  This 
test  showed  all  eight  saaples  were  randoa  at  the  95* 
confidence  level. 


TABLE  IV 

Suaaary  of  ABOVA  Basalts  at  95%  Confidence 


I 

I 


Rounded 

Data 

Unrounded  Data 

F* 

j 

nail 

1  null 

F  1 

hyp. 

F  |  hyp. 

300  s  arc  I 


isstrunent  j 

3.97  | 

O 

e 

o 

accept 

0.09 

accept 

observer 

2.70  ! 

| 

1.69 

accept 

1.69 

accept 

| 

interaction  | 

2.  70  | 

1.36  | 

accept  i 

1. 35 

accept 

50  0  a  arc 

instruaent 

3.92 

0.01  j 

accept 

0.01  I 

accept 

observer 

2.68 

0.13  j 

accept 

0.  14  | 

accept 

interaction 

2.68 

0.25  | 

accept  j 

0.25  | 

1 

accept 

700  a  arc 

instruaent  j 

|  3.  84 

l 

1  0.05 

1  accept 

0.05  I 

accept 

observer  1 

2.60  j 

j  1  .86 

i  accept 

1.86  j 

accept 

interaction  j 

2.60  j 

|  0.28  i 

1  accept 

0.  29  | 

accept 

1000  a  arc 

instruaent  | 

3.92  1 

0.01  | 

accept  | 

j 

0.  01 

accept 

observer 

2.68  ] 

0.45 

1 

accept  i 

0.45 

accept 

interaction  ; 

2.68  | 

1  .74  j 

accept  i 

1.79 

accept 

1500  a  arc 

instruaent 

3.  84 

0.78  i 
< 

l  accept  | 

0.  73 

accept 

observer 

2.68 

1  .44  j 

accept  | 

1.  41 

accept 

interaction 

2.68 

5.01  | 

i  reject  1 

5.02 

reject 

3000  ■  arc  I 


instruaent 

3.95 

0.10 

accept 

0.  12  I  accept 

observer 

2.70 

1  .69 

accept 

1.74  j  accept 

interaction 

2.70 

1.36 

accept 

1.37  accept 

C.  IHTEBPOLATIOH  BTALOATIOS 


The  two  methods  of  interpolation  used  by  N3S  were  next 
evaluated.  As  explained  in  section  III.C. ,  three  sets  of 
positions  were  computed.  These  were  a  set  of  actually 
observed  positions  using  Aztrac,  and  two  sets  of  corre¬ 
sponding  positions  computed  by  the  two  different  interpola¬ 
tion  algorithms  discussed  in  section  II. C.  These  two 
algorithms  were  evaluated  by  computing  the  distance  between 
each  observed  position  and  each  corresponding  interpolated 
position.  A  PORTBAH  program  written  by  the  author  performed 
the  computations.  The  NGS  geodetic  direct  subroutine 
[Pfeifer,  1975]  was  used  to  compute  positions  from  observed 
distances  and  Aztrac  directions  to  the  vessel.  The  interpo¬ 
lated  positions  were  computed  using  equations  2.23  and  2.24, 
which  are  the  same  algorithms  used  in  the  NOS  interpolation 
subroutine  TCARC  [Ehrhardt,  1979  ].  Distances  between  the 
two  corresponding  positions  were  computed  on  a  plans,  rather 
than  using  a  geodetic  computation.  There  is  negligible 
difference  between  plane  and  geodetic  methods  at  the 
distances  (about  ten  meters)  under  consideration  here. 
Results  of  these  computations  are  given  in  Table  V,  which 
shows  for  each  range  arc  the  average  distance  in  meters 
separating  the  interpolated  and  the  corresponding  observed 
positions.  An  indication  of  the  variability  of  these  values 
is  shown  by  the  percentage  of  interpolated  positions  falling 
farther  than  1.0  meter  away  from  the  actual  position. 

Little  can  be  inferred  from  the  results  in  Table  v 
because  this  table  is  really  an  intermediate  step  towards  a 
rigorous  evaluation  of  the  raw  data.  This  table  should  be 
viewed  as  only  a  general  indication  of  the  effectiveness  of 
interpolation,  since  it  is  a  direst  result  of  the  ability 
of  the  boat  operator  to  steer  the  vessel  in  an  arc,  wind  and 
sea  conditions  have  a  tremendous  effect  on  interpolation 
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TABLE  V 

Basalts  of  Interpolation  Evaluation 


Bange 

Angle 

Only 

Angle  and 

Distance 

a 

T 

*  >  1.0  m 
from  mean 

T 

*  >  1.0  m 
from  mean 

3000 

1.15 

47* 

2.41 

55* 

1500 

2.09 

85* 

4.98 

91* 

1000 

1.63 

72* 

4.07 

78* 

700 

2.04 

75* 

3.01 

80* 

500 

1.86 

79* 

4.  15 

33* 

300 

2.59 

79* 

5.32 

34* 

j 


effectiveness.  Although  this  experiment  was  carried  out  in 
representative  survey  conditions.  Table  v  should  not  be 
viewed  as  being  applicable  to  all  situations.  The  table 
does  indicate  that,  whenever  possible,  automatic  recording 
of  range  data  should  be  used. 

Pull  analysis  of  the  interpolation  algorithms  should  be 
the  two-dimensional  equivalent  of  tasting  for  the  difference 
between  means.  This  is  because  both  the  interpolated  and 
observed  positions  are  not  "true”  positions,  but  have  some 
error.  A  one-dimensional  test  of  differences  between  means 
is  well  established,  and  is  discussed  in  several  references, 
including  Honaacott  (1977)  .  The  null  hypothesis  for  such  a 
test  is 

d  -  (4.4) 
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where and ace  the  aeans  of  aha  two  populations,  and  d 
is  soae  arbitrary  distance  selected  by  the  experiaenter. 

The  two-diaensional  problea  has  tha  saae  null  hypothesis  but 
the  aatheaatics  of  the  test  have  not  been  established.  This 
problea  is  illustrated  in  Figures  <4.1  and  4.2. 


Figure  4.1  one-Dinensional  Difference  Between  Beans. 


A  proper  analysis  of  the  data  would  inquire  for  each 
interpolated-obsarved  pair  of  positions,  whether  the 
distance  between  the  two  positions  was  greater  than  d  for  a 
particular  confidence.  More  work  than  could  be  incorporated 
into  this  thesis  is  required  to  fully  evaluate  the  data. 


D.  AEALISIS  OF  FACIOHS  1FFECTIHG  THE  BBSOLTS 

The  results  of  this  experiment  ware  presented  in  Tables 
III,  IV,  and  V.  in  attenpt  will  now  be  aade  to  analyze  the 
experiment  for  errors  in  logic  and  technique,  in  order  to 
better  understand  these  findings. 


interpolated 


observed 


i 

_ j 


Figure  4.2  Two-Dimensional  Difference  Between  Means. 

The  *ost  important  results  from  this  thesis  are  the 
estimates  of  <r  for  theodolite  pointing  error.  It  is  obvious 
from  Table  III  that  standard  deviations  of  both  the  rounded 
and  unrounded  T-2  error  values  are  about  one-half  that  of 
the  Aztrac,  for  all  angular  speeds  considered.  The  ANOVA 
technique,  however,  shows  no  statistical  difference  between 
the  instruments  at  the  95*  confidence  level.  An  analysis  of 
the  data  used  to  obtain  the  results  yields  a  potential 
explanation  for  this  apparent  contradiction. 

The  original  data  (pointing  errors  in  seconds  of  arc) 
were  made  the  subject  of  empirical  probability  density  plots 
using  the  subroutine  HISTS  [Robinson,  1974]  on  the  UPS 
computer.  These  plots  show  probability  density  versus 
error,  as  well  as  mean  and  standard  deviation.  An  example 
of  these  plots,  for  the  500  meter  range  arc,  is  given  in 
Figure  4.3.  The  remaining  plots  are  found  in  Appendix  c. 
Plots  are  shown  for  both  Aztrac  and  T-2  (unrounded) ,  for 
each  range  arc.  A  striking  feature  of  the  Aztrac  curves  is 


1 - - - - - - - J 

Figure  4.3  Probability  Density  vs.  Error. 

their  bimodal  shape,  as  compared  to  a  single  peak  for  the 
T-2  curves.  It  can  be  easily  seen  from  the  curves  how  the 
spread  of  this  bimodal  distribution  would  increase  the 
computed  standard  deviation  for  the  Aztrac  data. 

The  method  of  data  acquistion  for  this  thesis  was 
seaiautcaated  in  that  all  T-2  angles  were  manually  recorded, 
while  the  Aztrac  angles  ware  recorded  by  pressing  a  button 
aboard  the  vessel.  It  is  probable  that  a  time  lag  existed 
between  all  the  T-2  observations  and  the  Aztrac  observa¬ 
tions,  despite  the  best  efforts  of  the  observers,  because 
the  observation  procedure  was  not  totally  automated.  If 
this  were  true,  there  would  be  little  difference  between  the 
observed  angle  e  for  the  T-2  and  the  computed  angle  0C  , 
because  ec  is  associated  with  a  reference  position  also 
derived  froa  T-2  observations.  The  observed  angle  e„  for 
the  Aztrac  would,  however,  be  consistently  different  from  9e 
because  of  this  time  lag  and  because  the  vessel  was  moving 
to  the  left  or  right  with  respect  to  the  Aztrac  observer. 
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Since  approximately  equal  numbers  of  observations  were  made 
with  the  boat  moving  to  the  left  or  right  along  the  range 
arc,  the  distribution  of  Aztrac  pointing  errors  would  take 
on  a  bimodal  3hape. 

It  is  understood  that  the  data  analyzed  by  these  curves 
come  from  not  one  but  four  different  samples,  and  than  the 
curves  should  not  be  expected  to  be  perfectly  peaked.  The 
T-2  curves,  however,  also  come  from  four  samples  and  do  not 
have  multiple  peaks.  This  analysis  is  further  supported  by 
finding  the  distance  between  one  peak  of  the  Aztrac  curve 
and  the  single  peak  of  the  T-2  curve  in  figure  4.3.  The 
distance  in  arc  seconds,  when  converted  to  meters,  is 
roughly  the  distance  the  vessel  traveled  in  one  second.  One 
second  of  time  is  certainly  a  reasonable  figure  for  the  time 
lag  discussed  above.  A  manual  check  of  the  raw  data 
recorded  in  the  field  also  suggests  such  a  time  lag.  The 
original  data  were  sorted  into  two  sets  of  ‘•left**  and 
"right"  observations,  which  were  analyzed  for  mean  and  stan¬ 
dard  deviation.  Results  of  the  analysis  are  shown  in  Table 
VI.  This  table  gives  the  mean  and  standard  deviation  for 
the  "left"  and  "right"  data  sets,  and  shows  chat  "the  mean  of 
both  sets  was  about  two  meters  to  the  left  or  right  of  the 
reference  position  of  the  vessel.  This  two  meter  difference 
corresponds  closely  to  a  nominal  vessel  speed  of  two  meters 
per  second  (four  knots)  for  the  boat  used,  and  a  time  lag  of 
one  second,  deans  for  the  1500  and  3000  merer  ranges  are 
somewhat  unequal  because  sea  conditions  at  these  offshore 
ranges  caused  the  boat  to  travel  slower  in  one  direction. 

The  rightmost  column  in  Table  VI  gives  the  pooled  standard 
deviation  of  each  "left"  and  "right"  data  set,  which  is  the 
best  estimate  of  the  population  standard  deviation  ,  <Ta  ,  for 
the  Aztrac. 
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TABLE  71 

Corrected  Aztrac  Standard  Deviation 

{  {  Kean  ,  T  1  s  I Aztrac  Pooled 

Arc  jleft/J  i  Stajjda^d 

(right  i  a  (sec)  |  a  (sec)  |  Deviation 

)  L 

300 

1  8 

2.95  (  20  35) 
-2.  67  (-1836)  | 

1.17  (  806) 

0.99  (682) 

1.07  (741) 

500  i  B 

2.  40  (  993)  | 

-2. 69  (-11  12) 

1.50  (660) 
1.31  (539) 

1.44  (598) 

700  1  B 

1.61  (  476) 

-2.37  (  -700) 

0.90  (266) 
1.28  (379) 

1.05  (312) 

1000  i  J 

2.49  (  514) 

-2.92  (-6  04) 

1.55  (321) 

1.30  (268) 

1.43  (296) 

1500  1  L 

1  8 

1.93  (26  6) 

-3.  22  (-444)  i 

1. 66  (228) 
1.06  (147) 

I  1.45  (199) 

3000  j  J 

0.48  (33)  1 

-3.  97  (-273) 

1.42  (98) 
1.22  (84) 

I  1.30  (90) 

The  systeaatic  error  caused  by  a  tine  lag  as  the  boat 
aoved  left  or  right  in  the  observer's  field  of  view  was  the 
result  of  faulty  design  of  the  experiaant.  Tha  proper  way 
to  correct  this  problea  would  be  to  duplicate  the  experiment 
using  better  synchronization  of  all  observations,  in  alter¬ 
native  would  be  to  aodel  the  systeaatic  error  and  apply 
corrections  to  the  existing  data.  Such  a  aodel  should 
include  an  estiaate  cf  the  boat  speed  and  its  left  or  right 
direction  with  respect  to  the  observer. 

Different  AN07A  results  night  be  obtained  using  the  data 
corrected  for  systeaatic  error,  but  this  would  not  explain 
the  A 9071  results  in  Table  17.  Tk.  conditions  affecting 
data  acquistion  east  again  be  considered,  as  well  as  an 
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understanding  of  the  ANOVA  process,  when  offering  an 
explanation. 

The  denominator  of  the  F-ratio  discussed  in  section  B  is 
essentially  the  experimental  error  of  the  measurement 
process.  Thi3  error  is  primarily  a  function  of  the  error 
<  ellipse  for  the  NreferenceH  position  shown  in  Figure  i.3. 

Table  VII  gives  an  estimation  of  this  experimental  error,  by 
comparing  the  size  of  the  error  in  B0  and  0C  .  Uncertainty 
in  6C  is  given  as  the  mean  major  axis  of  all  error  ellipses 
for  a  given  range  arc.  Uncertainty  in  0O  is  the  pooled 
standard  deviation  ,  ,  of  each  te3t  theodolite.  The 

values  for  in  the  Aztrac  column  are  from  Table  VI,  and  the 
T-2  values  are  from  Table  III.  An  examination  of  Figure  3.3 
shows  that  use  of  the  major  axis  of  the  ellipse  is  a  worst 
case  estimate  of  the  uncertainty,  since  the  ellipse  could 
have  any  orientation  in  the  x-y  plane.  Thus  it  can  be  seen 
from  Table  VH  that  the  uncertainty  in  6t  is  smaller  than 
that  of  the  observed  azimuth  6«  .  This  comparison  is  an 
indicator  of  the  precision  of  the  experiment. 

If  this  error  in  the  computed  angle,  0fc  ,  could  be 
reduced  by  decreasing  the  size  of  the  error  ellipse,  the 
denominator  of  the  P-ratio  would  be  smaller  and  the  ratio 
itself  would  be  larger.  Thus  a  more  precise  experiment 
could  produce  a  rejection  of  the  null  hypothesis,  although 
this  is  net  indicated  in  light  of  the  values  of  tfe  for  the 
Aztrac  and  T-2  shown  in  Table  VII.  Seducing  experimental 
error  any  further  than  this  study  would  be  difficult  under 
typical  hydrographic  conditions,  because  the  three  LOP 
theodolite  intersection  position,  adjusted  by  the  least- 
squares  method,  is  one  of  the  most  precise  positioning  means 
available  today.  The  size  of  the  reference  position  error 
ellipse  could  possibly  have  been  reduced  if  better  intersec¬ 
tion  angles  were  available,  but  this  was  not  possible  with 
the  particular  geographic  shape  of  aonterey  Bay. 
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TABLE  TII 

Experiment  Precision  and  Theodolite  Error 


Range, 

meters 

Mean  major  Aztrac  pooled 

axis  std. deviation 

aeters  (sec)  meters  (sec) 

<si  <r. 

T-2  pooled 
std.  dev. 
meters  (sec) 

tfx  <re 

300 

.974  (668) 

1.07  (741)  | 

1.26  (868) 

500 

.838  (344) 

1.44  (598)  | 

1.72  (712) 

700 

.754  (222) 

1.05  (312) 

1.30  (382) 

1000 

.824  (170) 

1.43  (296)  | 

1.08  (225) 

1500 

.738  (102) 

1.45  (199) 

0.9  2  (126) 

3000 

1.188  (  82) 

1.30  (90)  | 

1 

1.41  (97) 

In  sauary,  the  analysis  of  variance  performed  or.  these 
data  do,  in  feet,  shoe  that  there  is  no  statistical  differ¬ 
ence  between  iztrac  and  T-2  at  95X  confidence.  This  proce¬ 
dure  is  a  strictly  numerical  one  and  must  be  viewed  in  light 
of  the  original  research  conditions.  When  these  conditions 
are  carefully  analyzed  to  remove  as  much  systematic  error  as 
possible,  the  data  strongly  suggest  that  the  A  HD  Vi  results 
are  indeed  correct. 

E.  APPLXCATIOI  TO  FOSXTXOI  BBBOB  STAIDAB DS 

Prom  Chapter  I  it  can  be  seen  that  there  is  some  confu¬ 
sion  in  the  hydrographic  community  as  to  the  application  of 
probability  to  positioning  standards,  iith  this  in  mind, 
there  appear  to  be  four  possibilities  for  consideration  as 
standards  with  which  to  compare  tha  results  of  this  thesis, 
as  given  below. 

(i)  The  1968  IHO  standard  of  1.5  am  at  the 
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survey  scale,  with  the  90%  probability 
suggested  by  Hunson  (1977). 

(ii)  The  1982  IHO  standard  of  1.0  on  at  the 
scale  of  the  survey,  with  90%  probability. 

(iii)  The  current  range-azimuth  standard  of  NOS 
(0.5  an  at  the  survey  scale)  assuming  a 
probability  associated  with  d^mS . 

(iv)  The  drm5  standard  of  microwave  range- 
range  positioning  found  in  Umbach  (1976). 

This  requires  that  drw»i  values  at  the 
survey  scale  not  excead  0.5  an  for  1:20,000 
scale  surveys  and  smaller,  1.0  am  for 
1:10,000  scale  surveys,  and  1.5  mm  for 
surveys  of  1:  5,000  scale  and  larger. 

These  four  possibilities  nay  now  be  compared  to  the 
position  errors  of  the  Aztrac  and  1-2  by  using  the  pointing 
error  results  given  in  Table  VII.  Since  the  precision  of  a 
range-azimuth  positicn  is  also  a  function  of  the  standard 
deviation  of  ranging  error,  ,  assumed  values  of  1.0  and 
3.0  neters  are  used  in  this  analysis.  Some  manufacturers  of 
nicrovave  ranging  equipment  used  in  range-azimuth  posi¬ 
tioning  report  a  1.0  meter  value  for  ,  but  3.0  meters  is 
nost  often  used  by  NOS  personnel  [Wallace,  1983]  and  has 
sone  supporting  experimental  evidence  by  unbiased 
experimenters  [  Munson,  1  977  ]. 

Using  these  assumed  values  for  (T,  and  the  observed 
values  for  <Tt  given  in  Table  VII,  the  error  ellipse  axes  may 
be  computed  using  equations  2.11  and  2.12.  However,  the 
error  ellipse  must  be  converted  to  a  circular  error  figure 
in  order  to  be  compared  with  the  standards  above.  The 
assumed  standards  (i)  and  ( ii)  must  use  a  90%  probability 
circle,  which  can  be  computed  using  Table  1,  and  standards 
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angular  speed  (arc  min/ sec) 


Pigare  4.4  Basalts  coapared  to  d . 

(Hi)  and  (it)  aust  ass  ths  drmt  for sala  given  by  equation 
2.16.  Basalts  of  thsss  cos  potations  ara  shown  graphically 
in  Pigura  4.S  for  90*  probability,  and  in  Pigara  4.4  for 
drm, .  Each  figure  has  angalar  spaad  of  tha  thaodolita  along 
the  abscissa,  and  an  ordinate  consisting  of  a  distance 
scale,  in  aetars,  indicating  tha  radius  of  tha  error  circle. 


The  error  circle  for  Figure  4.5  is  associated  with  the  90% 
probability  level,  while  the  drmj  circle  for  Figure  4.4  is 
of  soeewhat  variable  probability,  as  discussed  in  Chapter 
II.  The  points  plotted  in  Figure  4.5  have  eccentricity 
values  ( 0*/c,  >  ranging  froa  0.3  to  0.6,  which  indicate  a 
probability  range  of  froa  about  65%  to  67.5%. 

These  figures  also  clearly  show  an  ioproveaent  in  the 
estiaate  of  <T$  for  the  Aztrac  as  a  resain  of  the  time  lag 
correction  discussed  earlier.  The  position  error  values  for 
the  uncorrected  Aztrac  error  is  shown  by  triangles  in  both 
figures,  while  position  error  computed  using  the  corrected 
(f0  values  are  shown  by  solid  dots.  For  both  Aztrac  cases, 
two  solid  linear  regression  lines  are  drawn.  Both  figures 
show  that  the  uncorrected  Aztrac  values  have  a  much  greater 
variability  than  the  corrected  values,  and  that  the 
corrected  values  are  almost  the  same  as  the  T-2  position 
error  values  indicated  by  open  dots  and  dashed  linear 
regression  lines.  If  the  corrected  Aztrac  values  are  taken 
to  be  the  best  estiaate  of  position  error  for  this  instru¬ 
ment,  then  a  relatively  constant  error  is  indicated  for  the 
entire  range  of  angular  speeds  considered  here.  This  is 
about  3.3  meters  d^j  and  4.6  meters  (90%  probability)  for 
both  instruments  using  a  0",  value  of  3.0  meters,  and  1.6 
meters  drnl(  and  2.5  meters  (90%)  for  a  O',  value  of  1.0 
meters.  Plots  are  not  shown  for  computations  using  a  range 
error  of  1.0  aetar. 

The  four  possible  assumptions  for  position  error  stan¬ 
dards  are  compared  to  Aztrac  and  T-2  position  errors  in 
Table  Till.  The  roaan  numerals  heading  the  columns  of  this 
table  refer  to  the  position  standards  associated  with  the 
same  numerals  at  the  beginning  of  this  section.  The  reader 
should  use  the  table  by  selecting  one  of  these  columns  and 
inspecting  it  froa  the  top  to  the  bottoa  of  the  table.  The 
first  row  of  Table  Fill  indicates  the  probability  associated 
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Figure  4.5  Basalts  compared  to  90*  Probability. 


with  each  positioning  standard.  Tha  second  row  lists  tha 
saxiaus  position  error  allowed  by  that  standard,  at  the 
scale  of  the  survey.  Bows  three  and  four  show  the  errors 
allowed  in  row  two,  when  converted  to  actual  distances  for 
two  representative  survey  scales  of  1:5,000  and  1:10,000. 
Bovs  five  and  six  show  the  radius  of  the  associated 


probability  circle  fcr  both  Aztrac  and  T-2,  for  ranging 
error  values  of  3.0  and  1.0  aeters,  respectively.  The  dual 
probability  percentages  in  coluans  (iii)  and  (iv)  indicate 
the  variable  probability  of  d^,  .  Dual  values  in  cclunn 

(iii)  for  aaxiaua  error  at  the  survey  scale  result  from  the 
NOS  standard  for  range-range  positioning.  The  remainder  of 
the  table  presents  conclusions  as  to  whether  the  T-2  and 
Aztrac  meet  the  various  standards.  For  example,  in  column 

(iv)  the  observed  3.3  aeter  drM,  value  in  row  five  is  less 
than  the  aaxiaua  allowable  error  of  5.0  aeters  shown  in  row 
four.  Therefore,  the  T-2  and  Aztrac  do  meet  th9  NOS  range- 
aziauth  standards  of  0.5  an  at  the  scale  of  the  survey  for 
1:  10,000  scale  surveys. 


Position  Standards  Comparison 
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Of  the  original  objectives  for  this  thesis  discussed  in 
Chapter  I,  the  first  and  most  basic  was  the  deter a inarion  of 
pointing  error  standard  deviation  for  the  Aztrac  and  T-2 
theodolites.  No  investigation  of  this  type  had  ever  been 
done  for  conditions  typical  of  a  range-azimut h  survey.  An 
experinent  was  carefully  designed  to  deternine  this  pointing 
error  and  to  deteraine  if  there  was  a  statistically  signifi¬ 
cant  difference  between  the  instruaents.  The  initial  esti- 
aates  of  pointing  error  were  given  in  Table  III,  which  shows 
the  Aztrac  to  have  an  error,  when  converted  to  distance,  cf 
about  3.0  aetars,  while  the  estiaate  was  about  1.3  neters 
for  the  T-2 . 

An  uncoapensated  systeaatic  error  in  the  data,  due  to 
the  tiae  lag  discussed  in  Chapter  III,  was  discovered  when 
eapirical  probability  density  function  plots  were  made  of 
the  entire  data  set  for  each  instrument.  This  led  to  a 
revised  estiaate  of  the  pointing  error  for  the  Aztrac, 
because  the  biaodal  distribution  caused  by  this  time  lag 
adversely  affected  the  original  estiaate  for  <f9 .  This 
revised  estiaate  is  about  1.3  meters,  as  shown  in  Table  VII, 
and  is  almost  the  saae  as  the  value  of  <f9  for  the  T-2 
instrument.  When  these  estimates  for  C9  are  viewed  in  light 
of  the  precision  of  the  experiment,  as  shown  in  the  table, 
it  is  seen  that  the  actual  values  of  could  be  smaller 
than  indicated,  because  smaller  values  would  be  masked  by 
the  relative  imprecision  of  the  experiment.  It  is  to  be 
concluded,  however,  that  the  actual  values  of  the  pointing 
error  of  each  instrument  are  no  larger  than  those  given 
here. 
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The  question  of  a  statistically  significant  difference 
between  the  instruments  was  then  considered  using  the  A NO VA 
technique,  which  can  be  said  to  coapara  a  variance  component 
due  to  the  instruments  with  a  variance  component  due  to  the 
precision  of  the  experiment.  This  precision  was  not  very 
much  greater  than  the  variance  of  the  instruments,  but  was 
based  on  the  most  precise  positioning  method  generally 
available  for  hydrography  --  an  intersection  position  using 
three  theodolites.  The  ANOVA  procedure  indicated  no  signif¬ 
icant  difference  between  instruments,  but  if  the  precision 
of  the  experiment  had  been  increased,  a  significant  differ¬ 
ence  could  possibly  have  been  detected.  In  light  of  the 
subsequent  discovery  and  elimination  of  the  systematic  error 
due  to  a  time  lag,  this  conclusion  of  no  difference  between 
theodolites  appears  to  be  well  justified. 

An  evaluation  of  interpolation  methods  was  the  second 
objective  of  this  thesis,  and  although  the  analysis  was  not 
as  rigorous  as  it  could  be,  it  can  be  concluded  that  there 
is  a  measurable  distance  between  an  interpolated  position 
and  a  corresponding  observed  position.  This  has  never  been 
done  for  the  case  of  range-azimuth  positioning  because  the 
rapid  position  fixing  available  with  Aztrac  has  not  been 
available.  It  has  been  shown,  through  an  error  propagation 
analysis  of  the  interpolation  algorithms,  that  the  interpo¬ 
lated  error  is  not  inherently  due  to  the  algorithms  them¬ 
selves.  The  arror  is  therefore  due  to  the  inability  of  the 
vessel  tc  follow  the  range  arc,  which  is  caused  primarily 
because  of  environmental  conditions  and  the  vessel  opera¬ 
tor*  s  track  keeping  capability.  The  distance  between  the 
interpolated  and  corresponding  observed  positions  may  be  as 
much  as  two  to  four  meters,  as  indicatad  in  Table  V,  and  is 
roughly  twice  as  great  for  a  position  that  is  computed  using 
both  distance  and  angular  interpolation,  as  for  a  position 
using  angle  interpolation  alone.  It  is  therefore 
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recommended  that,  whenever  possible,  automatic  recording  of 
range  data  should  be  used. 

The  third  and  most  important  objective  of  the  thesis  is 
a  comparison  of  the  total  position  arror  using  these  instru¬ 
ments  with  the  required  error  standards  of  the  major  hydro- 
graphic  survey  organizations.  The  lower  half  of  Table  VIII 
gives  these  results,  which  are  that  all  the  standards 
considered  are  indeed  met,  except  the  NOS  range-azimuth 
standard  at  1:5,  000  scale  using  a  range  error  of  3.0  meters. 
This  conclusion  requires  a  very  important  qualification 
regarding  the  T-2  instrument  and  the  errors  encountered 
while  pursuing  the  first  two  objectives.  These  are  errors 
due  to  interpolation,  and  to  the  time  lag  discussed  in 
Chapter  IV. 

The  approximately  one  second  time  lag  discovered  with 
the  Aztrac  data  set  is  not  actually  associated  with  the 
Aztrac  at  all  but  is  associated  with  the  T-2.  It  appeared 
to  be  a  systematic  error  of  the  Aztrac  in  this  experiment 
only  because  the  reference  postions  were  obtained  using  T-2 
instruments  similar  to  the  test  T-2.  It  must  be  concluded 
from  the  data  acquired  in  this  project  that  there  exists, 
for  any  angle  measured  with  a  T-2,  a  time  lag  of  about  one 
second  between  angle  observations  and  any  measurement  made 
aboard  the  vessel,  including  both  automatic  and  manually 
recorded  depth  and  range  data.  There  is  then  an  associated 
position  error  for  these  measurements,  the  magnitude  of 
which  depends  upon  vessel  speed,  which  was  about  two  meters 
for  the  four  knot  speed  used  in  this  experiment.  The 
conclusions  and  position  accuracies  for  observed  T-2  posi¬ 
tions  in  Table  VIII  do  not  take  this  additional  error  into 
account,  when  the  error  contributions  from  both  the  time 
lag  and  interpolation  are  considered,  it  can  be  concluded 
that,  positions  interpolated  between  observed  T-2  positions 
have  an  additional  error  of  about  two  to  four  meters.  Thus 
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the  total  actual  position  error  for  T-2  positions  and  posi¬ 
tions  interpolated  between  T-2  positions  sight  fail  to  meet 
more  of  the  standards  than  are  indicated  in  Table  VIII. 

Hawing  reached  conclusions  related  to  the  objectives  of 
the  thesis,  another  set  of  conclusions  and  recommendations 
can  be  made  for  the  Odom  iztrac  theodolite,  regarding  its 
ease  of  use  and  suitability  for  range-aziauth  hydrography. 
The  Aztrac  instrument  was  expressly  designed  for  range- 
azimuth  or  azimuth-azimuth  positioning,  and  has  features 
which  are  advantageous  to  the  operator  of  the  instrument  in 
the  field.  Such  advantages  are  discussed  in  Chapter  I,  and 
include  ease  of  tracking,  because  of  an  upright  telescope 
image  and  an  infinitely  geared  tangent  screw. 

The  Aztrac  theodclite  possesses  advantages  much  more 
important  than  ease  cf  use  in  the  field,  and  these  addi¬ 
tional  advantages  are  derived  primarily  from  its  ability  to 
be  interfaced  with  a  computer  aboard  the  survey  vessel. 
Besides  eliminating  systematic  error  due  to  a  time  lag,  a 
computer  based  survey  system  offers  the  additional  advantage 
of  being  able  to  measure  and  record  a  position  every  few 
seconds.  This  allows  three  important  advantages  over  a 
system  that  can  only  measure  positions  once  per  minute. 
First,  since  each  position  is  individually  measured,  no 
interpolation  i3  required  and  thus  better  accuracy  is 
obtained  than  with  a  nonautoaated  system.  Second,  no  manual 
data  logging  is  required,  which  reduces  blunders  and  greatly 
increases  the  speed  with  which  a  survey  may  be  processed. 
Third,  an  automated  system  allows  the  surveyor  to  run 
straight  sounding  lines  rather  than  curved  arcs,  because  an 
automated  system  can  provide  an  almost  real-time  cross  track 
error  indication  to  the  helmsman.  Sunning  straight  sounding 
lines  increases  survey  efficiency  by  orienting  the  lines 
sore  normal  to  the  depth  contours,  and  by  requiring  fewer 
total  linear  miles  of  hydrography  for  each  survey.  A  vessel 


nay  of  course  steer  straight  lines  with  a  nonautomated 
system,  but  it  is  extremely  difficult  to  maintain  the  strict 
line  spacing  requirements  for  hydrographic  surveys  without 
an  indication  of  cross-track  error,  so  curved  range  arcs  are 
usually  followed. 

Disadvantages  encountered  during  this  experiment  include 
additional  operator  fatigue  caused  by  the  requirement  to 
constantly  track  the  vessel,  '.’his  is  necessary  if  the 
substantially  increased  data  rate  available  with  this 
instrument  is  to  be  utilized.  The  extra  effort  to  track  the 
vessel  is  mere  than  offset,  however,  by  the  elimination  of 
manual  data  logging.  Care  is  required  by  the  operator  when 
rotating  the  instrument  through  an  arc,  because  if  the 
instrument  is  moved  too  rapidly  the  maximum  telemetry  data 
rate  is  exceeded  and  erroneous  angle  data  will  be  trans¬ 
mitted.  This  can  only  be  detected  by  checking  the  original 
initial  pointing.  Although  this  problem  was  observed  during 
a  manufacturer's  demonstration,  it  did  not  occur  during  the 
experimental  field  work.  Pinally,  the  transmitter  range  of 
5  km  is  rather  short  for  the  distances  used  by  NOS,  which 
can  be  up  to  about  10  km.  The  manufacturer  has  stated  that 
the  system  range  can  be  easily  extended  by  increasing  the 
transmitter  power  [Apsey,  1  983  ].  Although  this  thesis  meas¬ 
ured  Aztrac  error  at  a  maximum  distance  of  3  km  from  the 
shore  station,  the  conclusions  stated  here  should  not  be 
blindly  extrapolated  to  increased  ranges,  still,  if  the 
Aztrac  pointing  error  standard  deviation  is  reduced  to  its 
angular  resolution  (0.01  degree)  at  the  very  slow  angular 
speed  of  the  vessel  at  long  ranges,  the  error  appears  to 
remain  acceptable.  For  example,  0.01  degree  pointing  error 
at  a  range  of  10  km  results  in  an  error  of  only  1.74  meters. 

The  advantages  of  the  Aztrac  clearly  outweigh  its  disad¬ 
vantages,  so  it  is  therefore  recommended  that  the  Odom 
Aztrac  system  be  incorporated  into  the  computer  based  equip- 
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aert  used  by  NOS.  It  should  be  used  only  in  situations 
where  it  aeets  required  positioning  standards  based  on  the 
value  for  <f9  found  in  this  investigation.  It  is  also  recom- 
aended  that  the  systeaatic  error  inducad  by  the  tiae  lag 
discussed  here  should  be  accounted  for  in  operational  stan¬ 
dards  for  range-azimuth  positioning  using  the  siild  T-2. 

This  could  be  done  for  seai -automated  systeas  by  providing 
an  autoaatic  radio  signal  to  the  observer  on  shore  that 
precedes  the  actual  depth  measurement  by  one  second,  h 
simpler  method  of  reducing  this  systeaatic  error  could  be  a 
liait  on  the  speed  of  the  vessel,  depending  on  the  accuracy 
standard  required  for  the  particular  survey. 

A  final  recoaaendation  aust  be  aade  regarding  the  posi¬ 
tioning  standards  of  NOS.  At  present  there  are  conflicting 
standards  for  a  given  survey  scale,  depending  on  whether 
electronic,  hybrid  (including  range-aziauth)  ,  or  visual 
methods  are  used.  For  example,  positions  for  a  1:5,000 
scale  survey  aay  be  required  to  have  an  accuracy  of  either 
7.5  aeters  drwl  ,  or  2.5  meters  at  soae  unspecified  prob¬ 
ability,  depending  on  whether  microwave  range-range  or 
range-aziauth  methods  are  employed.  The  NOS  is  certainly 
the  most  progressive  hydrographic  organization  with  regard 
to  position  error  specifications,  but  it  is  recommended  that 
the  concept  of  probability  be  applied  to  all  positioning 
methods  and  not  only  electronic  ones.  Further,  if  the 
meaning  of  "seldom  exceed"  in  the  IHO  standards  is  to  be 
interpreted  as  a  90 %  probability  circle,  then  all  the  NOS 
standards  should  be  changed  to  reflect  that  standard. 
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BEABS  AND  STABDABD  DEflATIOIS  OP  ACQOIBED  DATA  SETS 


All  values  for  aeans  and  standard  deviations  are  in 
units  of  seconds  of  arc. 


TABLE  IX 
3000  Heter  Arc 
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TABLE  Z 
1500  Eater  Arc 
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T1BLS  XI 

1000  Heter  Arc 
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TABLE  XIII 
500  He  ter  Arc 
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TABLE  XIV 

300  Heter  Arc 
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GEODETIC  POSITIOH  OP  HOBIZOETIL  COSTHOL  STATIONS 


Station  Naae 

SOP  AR  (1947) 

OSS  HOH  (1978) 
90SSEL  (1932) 
AZTBAC 
T2 

GBOCEIVSB 


Lat  it  ida 


36 

36' 

32". 

,117 

36 

36* 

04". 

.685 

36 

37* 

18". 

,151 

36 

36* 

32". 

.530 

36 

36' 

32". 

,493 

36 

36* 

32". 

,512 

Longitude 


121 

53* 

24", 

.004 

121 

52* 

35" 

.9  00 

121 

54* 

11", 

.628 

121 

53* 

25" 

.310 

121 

53* 

25". 

.254 

121 

53* 

25". 

.286 

PCV) 
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BHPIHIC1L  PROBABILITY  DBHSIIY  PBICTIOI  PLOTS 


Figaro  Cl  Probability  Oansity  Plot:  300  a  arc. 


PC  V) 


Figure  C.2  Probability  Density  Plot: 


500  a  arc. 
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T-2:  SOLID  LINE 
flZTRfiC:  DASHED  LINE 


ERROR  IN  SECONDS 


Figure  C.  3  Probability  Density  Plot 
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Pigare  C. 4  Probability  Density  Plot:  1000  ■  arc. 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Or  STANDARDS  1%J  A 
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Figure  C.5  Probability  Density  Plot:  1500  ■  arc. 
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Figura  C.6  Probability  Oaosity  Plot:  3000  a  arc. 
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